NPS  ARCHIVE 
1968 
STEELE,  R. 


■ 


THE  THERMODYNAMIC  PROPERTIES  OF 
NITROGEN  -  OXYGEN  MIXmURFS 

AUTHOR:  LT  ROBERT  L.  STEELE,  USM 

THESIS  SUPERVISOR:  J.L.  SMImH,  JR. 
ASSOC.  PROF.  OF  MECHANICAL  ENGINEERING 

DATE:  17  MAY  1968 


DUDLEY  KNOX  LIBRARY 

NAVAL  POSTGRADUATE  SC 

MONTF.REY  CA  93943-5101  ^  ,„„ 

THE  THERMODYNAMIC  PROPERTIES  OF  "TE  SCHOOL 

.   NITRC  :   -OXYGEN  MIXTUR 
by 
ROBERT  LIVINGSTON  STEELE 
LIEUTENANT,  UNITED  STATES  NAVY 
B.S.,  United  States  Naval  Academy 
(I960) 
Submitted  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Naval  Engineer  and  the  Degree  of 
Master  of  Science  in 
Mechanical  Engineering 
at  the 
MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
May,  1968 


lUsh—£-&&h 


NAVAL  POSTGRADUATE  SCHOOL 
'  NTEREY," CALIF.  93940 

ABSTRACT 
Title:   The  Thermodynamic  Properties  of  Nitrogen-Oxygen  Mixtures 
Author:   LT  Robert  L.  Steele 

Submitted  on  May  17,  I96.8  to  the  Department  of  Mechanical 
Engineering  in  partial  fulfillment  of  the  requirements  for  the 
Master  of  Science  degree  in  Mechanical  Engineering  and  to  the  De- 
partment of  Naval  Architecture  and  Marine  Engineering  in  partial 
fulfillment  of  the  requirements  for  the  Professional  Degree, 
Naval  Engineer. 

This  thesis   i"nclii&e.s  a  survey  of  the  literature  through 
I967  on  the  methods  of  determination  of  the  thermodynamic  proper- 
ties of  nitrogen-oxygen  mixtures  in  the  two-phase,  liquid-vapor 
region.   The  virial  expansion,  empirical  equations  of  state,  pseu- 
docritical  constants,  excess  functions,  latent  heats,  and  graphi- 
cal methods  are  considered. 

The  methods  just  enumerated  are  evaluated  for  their  simpli- 
city, feasibility,  and  accuracy.   On  the  basis  of  the  state  of  the 
art  and  the* availability  of  appropriate  data,  thermodynamic  dia- 
grams are  constructed  for  1,  5»  10  a^d  20  atmospheres  o'f  pressure 
in  the  two -phase  region.   Diagram  coordinates  are  enthalpy- 
composition  and  entropy-composition.   The  diagrams  at  atmospheric 
pressure  show  good  correlation  with  known  data,  'while  those  at 
higher  pressures  are  best  estimates  using  excess  functions  and 
idealizing  assumptions. 
Thesis  Supervisor:   J.L.  Smith,  Jr. 
Title:   Associate  Professor  of  Mechanical  Engineering 
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LIST  OF  SYMBOLS 

a,  b;  am,  b^;  a^f  b^' single  substance,  mixture  and  pure  com- 
ponent of  mixture  constants  in  Redlich- 
Kwong  eqn.  of  state, 

a,  A,  b,  B,  ct  C  constants  in  Barker's  expression  for  G  . 

B(T),  C(T)  2nd  &  3rd  virial  coefficients  for  a  single  sub- 
stance . 

BiJt  B(T,x),  Bm 

---  2nd  &  3rd  virial  coefficients  for  a  mixture. 
C1Jk,  C(T,x),  Cm 

BAA»  BBB*  CAAA»  CBBB  2nd  &  ^rd  virial  coefficients  for  pure 

components  in  a  binary  mixture. 

BAB»  CABB»  CAAB  ^nd  &  3rd  interaction  or  cross  virial  coeffi- 

cients  for  a  binary  mixture. 

B(T),  C(T)  2nd  &  3rd  universal,  single  substance,  virial 

coefficients, 

b  "hard  sphere"  molecular  volume  =(2/3)7T  N0(j  ? 

Cg  specific  heat  at  constant  pressure  for  perfect  gas  or  at 

zero  pressure, 

D(T,x)  etc,  ^th  virial  coefficient, 

E  — --  internal  energy, 

f  fugacity  of  a  single  substance. 

f  ,  fugacity  of  component  i  in  solution 

f*  fugacity  of  pure  component  i  at  same  T,P  of  the  mixture. 

G  Gibbs  free  energy  =  H  -  TS. 

H  enthalpy  =  E  +  PV 

A  Hi  »  AK   differential  heat  of  vaporization  and  condensa- 

v  tion  at  constant  temperature. 

j  absolute  Joules 

k Boltzmann  constant  =  1.3805^  x  10~^  j/q^ 

L,  M,  N  defined  by  Eq.  (58). 

n  number • of  moles. 
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P  pressure. 

Pi  pure-component  pressure  of  component  i  in  a  gas  mixture, 

Pi  partial  pressure  of  component  i  in  ideal  gas  mixture, 

Pi  =  x\   P  (Dalton's  Law), 

R  universal  gas  constant  =  0.0820535  liter  atm/mole   K 

R  any  extensive  property 

r^ partial  property  of  component  i  in  solution. 

S entropy 

T  absolute  temperature, 

T   — -  base  or  reference  temperature, 

U  intermolecular  interaction  energy. 

V  molal  or  specific  volume 

V,  partial  molal  volume  of  component  i  in  solution, 

x^ mole  fraction  of  a  component  i  in  a  mixture  =  n./n. 

Z  ---  compressibility  factor  =  PV/RT 

superscripts 

E  -~-  excess  property, 

1  property  of  the  liquid  phase;  property  of  the  liquid  line, 

v  property  of  the  vapor  phase;  property  of  the  vapor  line, 

m  property  of  mixing. 

o  zero  pressure  state  or  perfect  gas  state  at  any  pressure, 

subscripts 

b  refers  to  the  normal  boiling  point. 

i  pure  component  property  of  a  mixture. 

m  property  of  a  mixture. 

c  critical  property. 

cm  pseudocritical  property. 


r  reduced  state, 

s  reference  state. 

Greek  symbols 

CX volatility  ratio. 

(j activity  coefficient. 

Q   general  function  of  temperature, 

£  characteristic  molecular  energy. 

Q   characteristic  molecular  distance, 

^ finite  change  of  quantity, 

[A — ■  -  chemical  potential  =  G 
All  state  properties  are  given  in  their  molal  or  specific  proper- 
ty form  unless  otherwise  noted. 
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INTRODUCTION 

It  is  the  purpose  of  this  thesis  to  present  the  results  of  a 
search  of  the  literature  and  the  subsequent  analysis  of  the  find- 
ings of  that  search.   The  subject  of  the  inquiry  was  the  thermo- 
dynamic properties  of  nitrogen-oxygen  mixtures  in  the  two-phase, 
liquid-vapor  region.   The  reasons  for  the  study  of  the  elemental 
properties  of  any  basic  substances  are  self-evident  especially 
for  substances  as  fundamentally  important  as  the  mixtures  of  nitro- 
gen and  oxygen.   It  should  suffice  to  say  that  the  hunt  carried 
the  author  through  the  fields  of  cryogenics,  low  temperature  phy- 
sics, high  pressure  physics,  thermodynamics,  physical  chemistry, 
molecular  chemistry  and  statistical  mechanics.   It  is  the  intent 
of  this  thesis  to  extract  from  these  areas  of  science  and  engi- 
neering the  most  appropriate  analytical,  graphical  and  machine 
methods  of  determination  now  available  and  to  present  the  proper- 
ties determined  therefrom. 

Acknowledgement  and  appreciation  are  expressed  for  the 

aid,  assistance,  and  direction  of  Professor  J.L.  Smith,  Jr. 

of  the  MIT  Cryogenics  Laboratory. 


PROCEDURE 

The  direction  of  this  work  consisted  of  two  distinct  phases. 
The  first  part  of  the  effort  was  directed  toward  a  thorough  search 
of  the  scientific  and  engineering  literature.   This  was  accom- 
plished by  a  personal  search  of  the  Chemical  Abstracts  through 
I967.   Machine  programmed  searches  were  conducted  for  the  author 
by  the  Cryogenics  Engineering  Laboratory  of  the  National  Bureau  of 
Standards  and  the  Defense  Documentation  Center.   These  initial 
efforts  were  principally  centered  on  nitrogen-oxygen  mixtures  and 
the  thermodynamic  approaches  to  the  derivation  of  properties  of 
binary  mixtures.  As  the  work  progressed,  molecular  theory  became 
of  interest,  and  the  extensive  bibliographies  in  the  major  texts 
of  the  field  were  utilized.   At  the  end  of  each  subsection  of  the 
Results,  the  references  are  listed  which  deal  explicitly  with 
that  subsection  as  annotated  in  the  text. 

The  second  part  of  the  program  was  the  analysis  of  the  data 
which  had  been  found.   It  must  be  said  that  the  determination  of 
the  properties  of  a  mixture  would  be  a  relatively  easy  problem  if 
the  mixture  behaved  simply  as  the  average  of  its  constituents. 
Unfortunately,  not  only  are  the  constituents  themselves  not  ideal, 
but  also  they  interact  with  each  other  in  a  non-ideal  manner. 

It  is  intended  in  this  thesis  to  discuss  the  PVT  behavior  of 
mixtures  and  the  various  methods  of  deriving  their  thermodynamic 
properties:   enthalpy  and  entropy.   Specifically  the  next  section 
will  cover  the  following: 

a,   the  theory  of  corresponding  states  and  pseudocritical  con- 
stants , 


b,  empirical  equations  of  state, 

c,  the  virial  expansion. 

d,  partial  property  correlations  and  mole  fraction  averaging, 

e,  excess  functions  as  deviations  from  ideality, 

f,  graphical  derivations, 

g,  latent  heat  techniques. 

No  attempt  has  been  made  to  make  interpretations  based  on  statis- 
tical mechanics  and  its  application  to  molecular  theory.   The  read- 
er interested  in  approaches  based  on  this  field  must  consult  the 
works  of  Eyring,  Kirkwood,  Lewinson,  Hirschfelder  et  al  for  their 
treatment  on  the  formal  and  model  approaches  in  this  area. 

The  final  section  of  the  Results  will  correlate  the  discussed 
data  and  provide  a  graphical  presentation  of  the  enthalpy  and  en- 
tropy of  binary  mixtures  of  nitrogen  and  oxygen  at  various  pres- 
sures. 


RESULTS 
Phase  Equilibria 

Perhaps  the  information  most  consistently  required  in  calcu- 
lations concerning  the  two-phase  region  is  the  equilibrium  behav- 
ior of  the  system  as  defined  by  Gibbs   1  ,   The  study  of  vapor- 
liquid  equilibrium  would  obviously  be  simplified  if  the  phases 
behaved  as  ideal  solutions,  for  then  we  should  only  require  infor- 
mation about  the  pure  components  in  order  to  develop  values  for 
the  mixtures. 

Deviations  from  ideal  behavior  occur  in  both  the  liquid  and 
the  vapor  phases  as  a  result  of  the  forces  that  exist  between  the 
molecules.   The  magnitude  of  the  intermolecular  forces  depends  on 
the  nature  of  the  molecules  and  the  average  distance  between  them 
as  expressed  by 

UAB(r)/€AB  =  f<r/CTab>  <X> 

This  expression  of  the  principle  of  corresponding  states  I  2,3J 
says  that  the  interaction  energy  of  a  pair  of  molecules  A  and  B 
at  a  distance  r  apart  may  be  expressed  in  terms  of  a  universal 
function  f  for  all  A  and  B,   The  parameters  £  AB  and  (j  AB  are  a 
characteristic  energy  and  distance  respectively  for  the  pair  of 
molecules.   In  the  liquid  state  the  molecules  are  relatively 
close  together,  and  the  force  fields  are  strong,  so  that  differ- 
ences in  shape  and  charge  distribution  will  cause  deviations  in 
behavior.   In  the  gas  phase  the  force  fields  are  not  as  strong, 
and  at  low  pressures  ideal  mixture  behavior  is  a  good  approxima- 
tion.  Mote  that  a  perfect  gas  is  defined  in  terms  of  negligible 


forces  between  molecules,  whereas  an  ideal  mixture  or  solution 
requires  only  the  same  forces  between  like  and  unlike  molecules. 
Hence,  at  low  pressures  ideal  solution  deviations  are  principally 
due  to  the  liquid  phase  while  at  high  pressures  both  vapor  and 
liquid  phases  contribute. 

It  is  perhaps  instructive  to  review  the  phase  behavior  of 
mixtures  under  pressure  [^,5,6J.   In  Fig.  1  the  curves  0CN2  and 
0Cq2  represent  the  usual  vapor  pressure  curves  for  the  pure  com- 
ponents terminating  at  their  respective  critical  points.   For  the 
binary  system,  however,  the  phase  boundary  is  the  region  enclosed 
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Figure  1.  Phase  boundaries  of  pure  components  and  mixtures 


by  the  curve  OC  0  within  which  both  the  liquid' and  the  vapor  coex- 
ist.  The  critical  point  of  the  mixture.  C  ,  is  defined  as  the 

m 

state  at  which  the  properties  of  the  vapor  and  liquid  phases  in 
equilibrium  are  identical.   If  the  phase  boundary  curves  of  sever- 
al mixtures  of  different  compositions  are  drawn,  the  locus  of 
their  critical  points  form  an  envelope  indicated  by  the  dashed 
line  in  Fig.  1.   The  maximum  pressure  and  temperature  for  the  ex- 
istence of  two  phases  are  called  the  cricondenbar ,  Cpf  and  the  cri- 
condentherm,  C™ ,  respectively   7  .   These  are  illustrated  in  Fig. 
2  as  is  the  phenomena  of  retrograde  condensation   8  . 
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In  Fig.  2  it  can  be  seen  that  if  a  line  of  constant  temperature 
is  followed,  condensation  will  begin  at  point  a  and  will  continue 
with  decreasing  pressure  until  maximum  condensation  is  reached  at 
point  b,  at  which  time  vaporization  will  occur  until  the  system  is 
again  all  vapor  at  point  c.   A  similar  phemomena  vrill  take  place 
at  constant  Pressure  with  vaporization  occurring  with  decreasing 


temperature  from  points  d  to  e,  • 

The  diagrams  most  frequently  used  in  representing  equilibri- 
um data  are  those  of  pressure  or  temperature  as  a  function  of 
molar  composition.   Fig,  3  shows  the  trend  of  the  phase  boundary 
curve  vrith  increasing  pressure  when  the  mixture  critical  tempera- 
ture and  pressure  fall  between  those  of  the  pure  components.   The 
area  of  two-phase  coexistence  decreases  with  increasing  pressure. 
The  phase  boundary  curve  breaks  away  from  the  pure  component  whose 
critical  temperature  is  lowest  and  recedes  along  the  critical  en- 
velope, finally  disappearing  altogether  at  the  higher  critical 
point. 

Turning  to  nitrogen-oxygen  mixtures  in  particular,  the  earli- 
est investigation  of  their  equilibria  was  conducted  by  Baly  in 
1900  at  atmospheric  pressure   9  #   The  temperature  scale  used  by 
him  was  in  error,  and  his  vapor  composition  measurements  proved 
to  be  innaccurate,   Inglis  I 10 J  in  1906  made  measurements  at 
7^.7  and  79.07°K  resulting  in  vapor  pressures  between  1  and  32 
atmospheres.   From  this  work  an  empirical  formulation  was  made 
for  the  tabulation  of  isotherms  and  isobars. 

During  the  next  few  decades  the  only  measurements  were  those 
of  Trapeznikova  and  Shubnikov  I  13J  in  193'+  on  an  isotherm  of 
85°K,  and  the  static  measurements  of  Sagenkahn  and  Fink  lM  at 
1  and  2  atm  in  19^. 

Little  else  was  done  until  1955  when  Armstrong  and  his  col- 
leagues at  the  National  Bureau  of  Standards  used  a  circulation 
method  to  make  accurate  determinations  of  the  compositions  and 
vapor  pressures  at  6$,    70,  and  77.5°K  [l5J .   Two  years  later 
Cockett  conducted  extensive  measurements  at  1,3158  atmospheres 


I l6jwlth  "results  reasonably  close  to  those  of  Dodge  and  Dunbar. 

Recent  significant  measurements  are  those  of  Din   17   v;ho 
used  the  static  method  for  nine  intermediate  mixtures  at  pressures 
to  10  atm.   His  results  likewise  show  those  of  Dodge  and  Dunbar 
to  be  reasonable  as  do  the  experiments  of  Latimer  1 18 ]  in  1957. 
In  196l  Yendall  and  Olzewski  followed  with  a  report  on  the  satur- 
ation properties  of  nitrogen-oxygen  mixtures  19J  . 

The  WADD  Technical  Report  of  1961  I  20 J  presents  a  compila- 
tion of  the  data  of  Armstrong  et  al,  Dodge  and  Dunbar  and  Sagen- 
kahn  and  Fink  on  T-x  diagrams.   The  data  of  Sagenkahn  and  Fink  at 
1  and  2  atm  and  of  Cockett  at  1.3158  atm  are  presented  on  P-x 
diagrams.   The  graphics  use  a  base  of  Hilsenrath's   21   pure  com- 
ponent vapor  pressure  values.   Tables  of  selected  experimental 
values  of  all  of  the  above  are  given. 

Air  Product  and  Chemicals,  Inc.  also  published  under  govern- 
ment contract  a  technical  report  I  22  with  a  P-T  diagram  of  vapor 
and  liquid  curves  based  on  the  data  of  Dodge  and  Dunbar  for  the 
two-phase  region  and  the  vapor  pressure  data  of  Hoge  for  oxygen 

23  and -Michel s  et  al,  Friedman  and  White,  Dodge  and  Davis,  and 
Henning  for  nitrogen  I 2^, 25 , 26 , 27J .   The  critical  envelope  is 
based  on  the  work  of  Kuenen  et  al  I  28, 29 J  and  can  be  considered 
precise  to  £2  K,   The  work  of  Kuenen  et  al  also  indicates  the 
presence  of  retrograde  condensation  over  small  ranges  near  the 
critical  points. 

The  most  recent  work  is  that  of  Haselden  and  Thorogood  in 
I963  at  high  oxygen  concentrations  I  30  which  showed  a  good  corre- 
lation with  the  earlier  works  of  Din  and  Latimer  below  96  ^mole 
oxygen.   They  demonstrated,  however,  that  the  relative  volatility 
falls  off  at  high  oxygen  concentrations.   This  is  alleged  to  be 


indicative  of  the  occurrence  of  association,  that  is,  the  weak 
chemical  interaction  between  the  oxygen  and  nitrogen  molecules. 

For  the  work  at  hand,  reliance  has  been  placed  primarily  in 
the  data  of  Dodge  and  Dunbar,  because  it  is  considered  to  be  con- 
sistent and  detailed  over  a  wide  range  of  values.   Their  data 
can  be  considered  accurate  to  to.l°K  at  pressures  near  atmospheric 
and  to  about  _1°K  at  pressures  above  seven  atmospheres   22, 
Saturation  values  for  the  pure  components  have  been  taken  for 
this  work  from  Stewart  for  oxygen   31  and  Strobridge  for  nitro- 
gen  32  ,  since  the  data  of  these  investigators  was  used  for  all 
other  pure  component  property  values,   Stewart1 s  empirical  vapor 
pressure  equation  was  fitted  to  the  data  of  Hoge  I  23 J  while  that 
of  Strobridge  fits  the  data  of  Friedman  and  White   25 1  and  Arm- 


strong 


33 


Individual  temperature-composition  diagrams  have  been  drawn 
for  1,  5 1  10  and  20  atm  and  are  presented  in  Appendix  A. 
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The  VirigQ  Equation 

The  virial  equation  of  state  is  simply  a  series  expansion  of 
the  compressibility  factor  in  terms  of  density. 

Z  =  PV  =  1  +  B(T,x)  +  C(T.x)    D(T,x)  +  ...  (2) 

RT  V         V%  T3 

Because  of  the  molecular  interactions,  the  virial  coefficients  are 
functions  of  both  temperature  and  composition  and  can  be  written 
as  [l] 

n   n 
B  =B(T,x)=£   £  x^B.  .      and  (3) 

n        n        n 
Cm  ■  C<T'X>  -T     E    ExixJxkClJ*  <*> 

i=A  j=A  k=A 
For  a  binary  mixture  these  become 


Bm  =  XABAA  +  -2xAxBBAB  +  xBBBB      and  (5) 

Cm  =  XA°AAA  +  5xAXBCAAB  +  3xAxBCBBA  +  XBCBBB  (6) 

where  x^  is  the  mole  fraction  of  component  A, 

Molecular  theory  tells  us  that  for  non-polar  molecules  (angle 
independent  molecules  possessing  no  permanent  charge  distribution) 
such  as  those  being  investigated,  calculation  of  their  virial  co- 
efficients depends  on  the  choice  of  proper  potential  functions. 
The  potential  function  represents  the  attractive  and  repulsive 
forces  between  molecules  and  many  are  found  in  the  literature 
I  2,  3  •   Although  these  forces  are  generally  exponential  in  char- 


TTT 


acter,  they  have  generally  been  simplified  to  reciprocal  power 
forms. 

The  most  popular  relationship  is  of  the  form 

U(r)  =  K  r~n-  K2r~m  (7) 


where  the  first  term  represents  the  repulsive  potential  energy  . 
and  the  second  that  of  the  attractive  contribution.   Where  n  and 
m  are  chosen  as  12  and  6  respectively  there  results  the  well  known 
Lennard-Jones  potential  (compare  with  Eq,  1), 


U(r)  =  ^£  ((j/r)12-(07r) 


(8) 


It  can  be  shown  k      that  by  using  the  Lennard-Jones  potential  the 
pure  component  virial  coefficients  can  be  expressed  as 


B(T)  =  f(TV,b  )   or   B*(T)  =  B(T)/b„=  f(T*) 
o 


o 


(9) 


where  b   is  the  "hard  sphere"  volume  defined  by  (2/3)7Tnl7^  and 


T  =  T/(6A),   Similarly 
C*(T)  =  C(T)/bn3  =  f(T*/CA) 


(10) 


Sqs.  (9)  &  (10)  have  both  been  tabulated  and  plotted. 

The  interaction  or  cross  virial  coefficients  are  signifi- 
cantly more  difficult  to  determine.  For  an  ideal  solution  obey- 
ing Amagat's  Rule   5 


baB=(baa+bbb)/2 


CAATH2CAAA+CRPm)/3    -    (BAA-BRR) 


'AAB" 


'BBB 


(11) 


CABB=  *  CAAA+ 2CEBB ) /3    -    ( -AA "BBB ) 
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which,  is  the  Lewis-Randall  rule  for  mixtures  of  simple  gases  and 
results  in  the  following  for  a  binary  mixture: 


Bx=xABAA+xBB3B     and 


(12) 


Cm=xACAAA+xBCBBB  "  3^AXB^BAA"BBB^ 

It  has  been  shown   6  I,  however,  that  as  B^ ( T ) /(jpj^    is  a 
universal  function  of  kT/£^A,  so  too  1S  ^AB^T'X^CJ'^B^ a  function 
of  ^T/^ab*   ^^e  TTieasuremen't  of  second  virials  in  mixtures  of  di- 
lute gases  where  the  third  virials  are  negligible  can  be  used  to 
verify  the  various  combining  rules  derived  on  the  basis  of  this 
relation.   The  most  popular  and  simple  of  these  are  the  Lorentz- 
Berthelot  rules   7 

^AB=(^AA^  BB>*    and    CrAB=«JAA+CrBB)/2  (13) 

which  is  based  on  London's  theory  of  dispersion  energy  |8,9tlOJ • 
These  rules  give  a  better  account  of  the  second  virial  coeffi- 
cient of  simple  gases  than  the  Lewis-Randall  rule,  yet  there  is 
real  doubt  k     whether  any  rule  will  yield  accurate  values  of  B.-o 
even  under  the  best  conditions  of  small,  non-polar  and  spherically 
symmetric  molecules.   To  make  matters  worse,  very  few  experimen- 
tal measurements  of  second  virial  coefficients  have  been  made  on 
mixtures  of  simple  molecules.   Newer  and  more  accurate  functions 
such  as  the  Kihara  potential  are  presently  coming  under  more  and 
more  study  [ll ,12, 13,1^1 . 

When  consideration  is  on  mixtures  other  than  dilute  gases, 
no  exact  relation  is  available  to  correlate  the  interaction  ener- 
gies.  Present  practice   6   is  to  extend  the  principle  by 


M- 


considering  an  average  interaction  energy 


<JWr)>  /C(x)  =  f  [r/CT(x) 


'AB 


(1*0 


While  the  virial  equation  provides  the  best  theoretical  basis  for 
estimating  the  PVT  properties  of  mixtures,  it  cannot  deal  satis- 
factorily with  any  but  low  pressure  gas  mixtures.   The  inability 
to  estimate  third  and  higher  virial  coefficients  generates  an 
overall  loss  of  accuracy. 

The  thermodynamic  properties  of  enthalpy  and  entropy  are  de- 
rivable from  the  virial  equation  I ^,15 i 16,1? J  .   From  the  First 
and  Second  Laws  for  a  single-component,  closed  system  dF=TdS  -  PdV 
and  H53+PV  whence  dH=TdS+VdP.   If  this  last  is  divided  by  dP  and 
T  is  held  constant,  the  relation  ( &  H/  d P)T=V-T(  d V/  &T) p  results. 
If  Z=PV/RT,  then  (  £>H/ 0 P)T= ( -RT2/P) ( b  Z/d   T)p  or  after  integration 

r   P2 
(H2-H1)T=-RT2/   (dZ/£  T)pd(ln  P)  (15) 


If  P1=0  and  Tr=T/Tc,  Pr=P/Pc ,  then 

P, 

I 
0 


(H-H°)/RTC=-T2/  (  ^Z/e)Tr)prd(lnPr)  (16) 


vrhere  (H-H°)  is  the  enthalpy  change  from  P=0  to  P=P  at  constant 
terraperature.   Then  for  mixtures   3  t  Fqs .  (2)  &  (15)  give 

(Hm-H°)/RT=l/v[Bm-TdBm/dT]  +  1/V2  [  Cm-T/2(dCm/dT)]  +  •••     (17) 

From  the  Maxwell  relation  (53/uP)T  =  (dv/dT)p  and  Z=PV/RT  we 
have  (  dS/ft  P)T=-R/P  [z+T(  b  Z/$  T)p]  .   Since  the  entropy  blows 
up  as  pressure  goes  to  zero,  it  is  convenient  to  use  the  artifice 
of  integrating  this  relation  as  a  real  gas  to  a  pressure  near 
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zero  where  properties  approximate  those  of  a  perfect  gas.   The 
equation  is  then  integrated  using  perfect  gas  properties  up  to  a 
convenient  reference  pressure,  P  ,  so  that  low  pressure  effects 


Then 


rp  fp 

(S-S_)/fc=/  (l-Z)d(lnP)+ln(P/P)-  T   (dz/dT)_  d(lnP)        (18) 
Jo  s  Jo  * 


on  the  entropy  mutually  cancel  l\ 

-P  rP 

(l-Z)d(lnP)+ln(P  /P)-  T/ 
'0  s     ^0 

or  in  reduced  terms 

fVT  rpr 

(S-SS)/RW  (1-Z)  d(lnPr)  +  ln(Prg/Pr)  -  T J   (  c)  Z/ 5  Tp)p^d(  lnPp ) 

(19) 

where  (S-S_)/R  is  the  entropy  change  from 

TfP0  to  T,P  or  T  .P-.  to  T  .P^.   Then,  Eq .  (2)  &  Eq.  (18)  give 

(Sm-S„J/R=-lnP  -T/v[dB  /dTl-  B2  /2V2+C_/2V2  -  T/2V2[dCTn /dT|   (20) 
in  ins        s     L^Jm  i_roj 

Values  can  be  determined  for  dB/dT  by  differentiating  curves  of 
B*(T)  as  a  function  of  T" ,  using  the  relation  dB  /dT  =(t  /k/bQ) 
dB/dT . 

0 

Finally  to  determine  H  and  S_  the  relations 

"  mm 

6 

H^ExjHJ"  and   Sms=  E  xiSis  -  R  E  lnxi  +  RlnPs  (21) 

i  i  i 

can  be  used,  where  H.  is  the  enthalpy  of  component  i_  at  P=0  and 
S  is  the  entropy  of  component  i_  at  reference  pressure,  P  . 
The  virial  equation  is  not  the  most  convenient  method  of  deter- 
mining the  thermodynamic  properties  of  mixtures  unless  computer 
facilities  are  available.  If  a  virial  equation  has  been  formu- 
lated for  PVT  properties,  however,  it  would  not  require  much 
additional  effort  to  program  the  derivatives  of  the  equation. 
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Pseudocritical  Properties 

It  was  stated  earlier  that  the  molecular  energies  of  inter- 
action could  be  described  by  their  dependencies  on  a  universal 
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function  as  expressed  by  Eq.  (1*0.   The  molecular  theories  of 


Prigogine  et  al 


W 


in  which  explicit  relations  for  the  parame- 


ters C  (x)  and  Q   (x)  are  developed  in  terms  of  their  pure  and 

mixture  forms,  ie,  ^  AA'  ^ab'  ^BB'   AA'   AB*   BB'  throuSh  the 
Lennard-Jones  or  other  m-n  potential  forms,  are  formally  equival- 
ent to  the  pseudocritical  methods  of  Kay  et  al  I  2J.   The  relations 
kT/£(x)  and  V/(j*(x)   are  equivalent  to  universal  functions  which 

may  be  replaced  by  T/T  (x)  and  V/V  (x) ,  where  T  (x)  and  V  (x) 

c  c  c         c 

are  pseudocritical  properties  of  the  mixture  I  3 J.   The  pseudo- 
critical  concept  is  based  on  the  existence  of  a  mythical  pure 
fluid  at  the  termperature  and  pressure  of  the  mixture  which  exhi- 
bits the  same  properties  as  the  mixture.   The  critical  properties 
of  this  fluid  are  the  pseudocritical  properties  of  the  mixtures. 

A  sizable  amount  of  literature  has  been  written  on  the  deter- 
mination and  use  of  the  pseudocritical  constants   2,^,5f6,7 I. 
Kay's  simple  mole  fraction  averaging  rules  I  2j  are  the  most  fre- 
quently used;  in  fact,  the  various  rules  for  pseudocritical  temp- 
erature all  reduce  under  various  simplifications  to  that  of  Kay, 

If  T   ET  (x),  Kay's  rules  are  given  by 
cm   c 


2_^  x.T  .    (22)    and   P     £_,  x. 


T  =  7,x.T  .    (22)    and   P     7  .  x  P  (23) 


j~  i  oi  cm    ^   l  ci 


For  T   i  Kay's  rule  produces  values  within  2%   of  other  rules  for 
0.5<TCA/T  <  2.0  and  0.5^  PCA/?CB  ^  2  [s].   For  values  of  Pc  dif- 
fering by  20;^  or  more,  the  modified  Prausnitz  and  Gunn  rule   9 
should  be  used. 

P  =Z   RT   /V  =R  )  ,  x .  Z„  .     )  ,  x .  T  .  /  )  ,  x  V  (2*0 

cm  cm  cm  cm   'V  i  cl  t_j  i  ci  (—»      i  ci 
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Kay!s  rule  has  been  verified  to  predict  the  true  critical  temp- 
erature of  N9-0?  mixtures  to  within  10  K. 

Using  the  pseudocritical  rules  for  mixtures  of  nitrogen  and 
oxygen,  the  author  computed  enthalpy  deviations  by  the  method  sug- 
gested by  Reid  &  Sherwood.   A  summary  of  this  method  f ollows : 

1,  Pseudocritical  constants  for  the  mixture  were  calculated  from 
the  f ollowing  pure  component  values  : 

Nitrogen  Oxygen 

T  126.26°K  [ll.l^]  15^.77°K  [l2,13] 

P  33.5^  atm  [  11, 1^1  50. 1**  atm  [12,13] 

V  90.07  cm3/mole  [15]  .075  cm3/mole  \  12*1 

Zc  .291  [8]  .292  [s] 

2,  Base  enthalpies,  H°  ,  for  nitrogen  and  oxygen  were  selected 
based  on  the  enthalpies  of  the  perfect  gas  at  1  atm  at  the 
normal  boiling  point  and  triple  point  respectively  11,12  . 

3,  The  ideal  state  enthalpies  at  the  various  isotherms  of  inter- 
est, H  ,  are  then  calculated  from 

J-  • 

H?  =  $o{\°   dT  '  (25) 

where  C°  =29.121  j/mole°K  and  C   =  29.359  j/raole°K  [ l6J . 

b.     At  a  particular  temperature  the  perfect  gas  mixture  enthal- 

pies,  H  ,  may  be  calculated  from 
m 

1 

5.   The  mixture  enthalpy  at  the  system  pressure  is  then 


where  the  isothermal  enthalpy  deviation  H   -  H°  can  be  deter- 
mined from  suitable  tables  I  17  . 
Reid  and  Sherwood  predict  errors  of  10$  in  the  vapor  phase  and 
20^  in  the  liquid  for  this  method  18  J.   This  author  carried  out 
calculations  at  8^-°K  and  1  atm  and  at  112°K  and  10  atra.   In  the. 
vapor  phase  for  the  case  of  8^°K  and  1  atm,  although  the  calcula- 
ted values  differed  from  those  obtained  using  a  mole  fraction 
average  of  the  pure  component  enthalpies  by  as  much  as  35  j/mole, 
the  results  were  linear  and  had  the  same  slope  (Fig.  *4-a),   in  the 
liquid  region  at  this  temperature  the  results  showed  a  small  neg- 
ative deviation  from  linearity  (experimental  evidence  indicates  a 
positive  excess  heat)  with  the  absolute  difference  increasing 
with  increasing  ^mole  02  (Fig.  ^b).   At  the  higher  pressure  the 
calculated  vapor  isotherm  was  no  longer  linear,  but  exhibited  a 
sizeable  negative  excess  heat.   Likex^ise  the  liquid  isotherm  at 
this  pressure  displays  a  somewhat  larger  contrary  behavior  than 
before.   All  isotherms  are  extrapolated  into  the  wet  vapor  region, 
These  calculations  seemed  to  be  undesirable  for  their  deviations 
from  known  behavior,  ie,  negative  excess  heat.   The  tables  offer, 
however,  a  relatively  simple  computational  technique  and  are  ap- 
plicable to  the  liquid  phase  as  well  as  the  vapor.   Perhaps  better 
values  of  the  pseudocritical  properties  would  produce  more  logi- 
cal results,  but  for  the  present,  simple  mole  fraction  averaging 
seems  to  provide  more  satisfactory  answers. 

While  tables  such  as  those  of  Lyderson  et  al  do  not  exist  for 

r 

entropy  deviations  from  the  perfect  gas  base,  a  relationship  can 
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Figure   4a. 
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Figure  4b. 
Figure  4a  &  4b.  Liquid  and  vapor  lines  determined  by  various 
methods.  Dashed  lines  join  pure  component  values". 
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be  written  between  properties  of  other  values  which  are  tabulated. 
Sq.  (IS)  can  be  written  18 J 

(S  -  Sg)/R=(H  -  H°)/RT  -  ln(f/Pg)  (28) 

when  P=P  and  values  of  both  (H  -  H°)  and  (f/Ps)  have  been  tabu- 
lated. 
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Empirical  Equations  of  State 


Perhaps  the  most  underdeveloped  of  the  techniques  for  the 
determination  of  the  properties  of  mixtures  is  that  of  a  workable 
equation  of  state,  particularly  in  the  two-phase  region. 

The  Benedict-Webb-Rubin  equation  ll  with  eight  constants 
has  been  fruitfully  tested,  but  unfortunately  is  applicable  only 
to  light  hydrocarbon  mixtures. 

Some  effort  has  been  made  to  determine  coefficients  for  the 
two  constant  Redlich-Xwong  equation   2   as  follows: 


P  +  a/[r0,5v(V+B)J 


Rr 


where  for  a  mixture   a 


(V  -  b) 

*     y     • 


(29) 


and   b 


m    L — »  3 


ar.d  also  a*  ---o.'l-278/T  ?'5P  .   and  b*  =  0. 0%7/T  .  P., 
i  rici        i  ij.  ^j- 

It  is  reported  that  some  success  has  been  obtained  in  using 
reduced  temperatures  and  pressures  derived  from  Kay's  rules  in 
the  generalized  Beattie-Bridgman  equation  [3  . 


Pr  =    (Tr/Vri)(1   "    °'°^T^vri) 


V    .    +   0.1867  [(1   -  0.03833)/V   . 
ri  l  riJ 

-(0.^758/Vr2)    (i  -  0.112?)/Vrl]  (30) 

Recent  efforts  |^»5|  have  been  fairly  successful  in  treating 
gas  mixtures  of  CO  and  N2  with  mixture  rules  for  the  Martin-Hou 


equation  of  state, 
5 


=  £ 


1=1 


f^v  -  b) 


-i 


whe 


re  f   B  At  +  BiT  +  Ciexp(k»T/Tc ) 


(3D 
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None  of  the  equations  of  state  have  been  extensively  tested  with 
mixtures,  with  the  possible  exception  of  the  BWR  equation  for 
light  hydrocarbons,  and  all  require  machine  handling.   In  the 
derivation  of  thermodynamic  properties  using  "qs,  (15)  &  (18) 
the  results  will  be  only  as  good  as  the  equation  used,  and  some- 
times worse,  because  of  the  errors  involved  in  the  determination 
of  the  derivatives  of  the  equation. 
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Partial  Properties 

It  was  mentioned  earlier  that  use  might  be  made  in  some  cases 
of  partial  properties  in  determining  mixture  properties.   Partial 
molar  properties  define  the  properties  of  components  in  a  mixture 
as  opposed  to  those  of  the  pure  components  and  represent  a  part 
of  the  total  property  of  the  particular  mixture.   In  particular, 
a  property  is  defined  as 

Ri  *  ^^/Oni)n^ifpfT    or        R  =  )_>n1Ri  (32) 

where  R  is  any  extensive  property  and  n.  is  the  number  of  moles 
of  component  i  in  the  mixture.   These  relations  are  governed  by 


2'P 


the   Gibbs-Duhem   equation 

Z_,  n      dJL    -    0  (33) 

i 

If  R  s  R/n  =  /_j  x^R.  ,  then  it  can  be  shown  |  1   that 

R A  =  R  -  xB(c)R/c)xB)   and   RR  =  R  -  \(^/()xA)  (3*0 

The  Lewis  relation   2   for  fugacity  as  it  is  related  to  the 
chemical  potential  can  be  expressed  for  a  mixture  as 

Mi  =  (  c)G/^)xi)x  >PjT  =  RT  lnf1  +  0i(T)  (35) 

where  (y  *(T)    is  some  function  of  temperature  only.   As  the  pres- 
sure approaches  zero,  f  approaches  P  as  opposed  to  the  chemical 
potential  which  blows  up  as  can  be  seen  from  (  d  lii^/  c^t'P)^   =   V*. 
At  constant  temperature  and  composition  Eq.  (35)  can  be  written 

M 

RT  lnff^/f  x^)  =  f  (Vi  -  V1)dP  s  RT  In  ^  (36) 

J  o 

where  (j  ^    is  defined  as  the  activity  coefficient.   For  ideal  sol- 
utions this  reduces  to  the  Lewis  &  Randall  rule  [2  J 

?  =  x*f   such  that  J1   =   1 .  (37) 

While  in  an  ideal  solution  it  is  only  necessary  that  the  inter- 
molecular  forces  be  the  same  for  all  components,  in  a  perfect  gas 
the  molecular  forces  must  be  zero.   In  this  case  the  pure  compon- 
ent fugacities  are  equal  to  the  pressure  and  Raoult's  law  results 
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f .  =  P.  =  xvP  (38) 

iii 

It  should  he  mentioned  that  the  Lewis-Randall  rule  generally 
provides  a  good  approximation  for  most  real  non-polar  gas  mix- 
tures at  any  pressure  as  long  as  the  component  in  question  is  sim- 
ilar to  the  other  component  or  present  in  excess  of  the  other  com- 
ponent  5  • 

If  the  fugacity  of  the  pure  components  in  the  mixture  can  be 
determined,  the  enthalpy  deviation  representing  the  non-ideality 
of  the  mixture  can  be  determined  from  the  relation 

d[ln(f./x.P)]  /£>T   PfX  =  (Hj  -  H°)/RT2  •  (39) 

The  mixture  enthalpy  then,  follows  from  the  partial  property  and 
perfect  gas  base  state  relation 


K  =  Lx,(H,  -  H?)  +  2_,x,  H? 
m    ^   i   i    l     ^     i   i 


(^0) 


where  the  pure  component  perfect  gas  values  are  determined  as  be- 
fore. 
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Excess  Functions 

Before  discussing  deviations  from  ideality,  it  would  do  well 
to  define  somewhat  more  fully  the  ideal  mixture.   One  definition 

ll  states  that  it  is  a  mixture  whose  components  have  chemical 
potentials  given  by 


jU  (P,T,x)  =  ^(P.T)  +  RT  lnxi 


(M) 


where  U   . (P,T)  is  the  pure  component  potential  at  the  pressure 
and  temperature  of  the  mixture.   It  follows  that  for  an  ideal  mix- 
ture 


x  lnx 


H 


m   l_,  x   S   *•  R  /_,  x  lnx 

i   i  i      i   i    x 

=  Lx4y  ' 


S  (^2) 


The  log  sums  are  referred  to  as  the  ideal  free  energy  of  mixing, 
Gm,  and  the  ideal  entropy  of  mixing  (Gibbs  paradox),  Sm.   The 
ideal  free  energy  of  mixing  is  always  negative  and  that  of  entro- 
py always  positive.   The  ide'al  enthalpy  and  molar  volume  of  mix- 
ing are  equal  to  zero. 

Excess  properties  may  be  defined  in  these  terms  as  the  dif- 
ferences between  the  values  of  the  properties  in  the  real  mixture 
and  the  values  of  the  properties  given  by  Eqs ,  (^2)  for  an  ideal 
mixture  of  the  same  composition. 
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It  has  been  seen  that  the  chemical  potential  may  be  related 
to  the  activity  coefficient,  so  that  -for  a  real  mixture  at  low 
pressures 

^.(P.T.x)  -/^(P.T)  =  RT  Inx^  Ji  (^3) 

The  partial  molar  properties  may  then  be  given  by 

H^  =  Hj[  -  RT2(dln/i;c)T)p 

B±   =  Si  -  R  lnxiJ'i  -  RT<c)  ln^/dT) 

V4  =  Vi  +  RT(5ln7i/^P)T 

The  excess  thermodynamic  properties  may  then  be  obtained  by  sub- 
tracting* the  ideal  quantities  from  the  partial  molar  quantities. 


=  rt  [xim/i 


HE   =   RT2    l^  x  A  din  J  1/5T)p 

SE  =   -R    L,  ^lnfx^j^)    -RT   L  x^  dlr\J±/  d^) 

VS   =   RT    Lxi(blnJ1/d?)T 


>   (^5) 


The  excess  functions  are  related  by 


GE  =  HE  -  TSE  =  HE  +  T(  dGE/e)T) 


(^6) 


Considerable  effort  has  been  made  in  the  past  decade  in  the  mea- 
surement of  excess  functions.   In  particular,  regarding  Np-CU 
mixtures,  Knobler  et  al  [2,3,^]  measured  HE  and  VE  at  77°K;  Wil- 
helm  &  Schneider  [5],  GE  at  77.  5°K;  Pool  et  al  [6.1,  G  at  83.8°K; 
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and  Duncan  &  Staveley   7  ,  G'J  at  63.1^°K.   In  addition,  values 
of  G   have  been  calculated  by  some  of  these  authors  from  the 
vapor  pressure  measurements  of  Dodge  &  Dunbar  I  8J  ,  Armstrong  et 
al  [9]  ,  and  Din  [lOj .   These  results  are  plotted  in  Fig.  32  for 
equimolar  liquid  mixtures. 

It  can  be  shown  111  that  the  variation  of  vapor  pressure  at 
constant  temperature  can  be  related  to  the  excess  Gibbs  free  en- 
ergy.  Vapor  pressure  and  composition  can  be  measured  in  an 
equilibrium  still,  a  large  number  of  which  have  been  devised 
[ll-2^J .   The  principal  problem  in  all  stills  is  the  difficulty 
in  eliminating  vapor  condensation  before  it  reaches  the  vapor 
sampling  point  and  of  preventing  liquid  droplets  from  being  car- 
ried over  to  the  same  point.   These  problems  can  be  eliminated  by 
methods  which  utilize  only  the  total  static  vapor  pressure  as  a 
function  of  liquid  composition.   Then  the  vapor  composition  may 
be  found  analytically  using  the  Gibbs-Duhem  equation.   Several 
methods  have  been  used  to  measure  the  total  vapor  pressure   25- 
29]. 

If  isothermal  vapor  pressure  measurements  are  known,  values 

E  1   v 

of  C"J  may  be  derived  from  knowledge  of  P,  x  ,  x  and  the  second 

virial  coefficient.   If  the  easier  experimental  procedure  has  been 

followed,  which  measures  x  and  not  xv,  there  are  several  graphic, 

algebraic,  and  machine  methods  which  may  be  used  to  solve  for  GJ 

[27,30-35J .   These  solutions  are  more  complicated  than  those  used 

when  both  compositions  have  been  measured.   Other  methods  have 

been  devised  for  the  calculation  of  the  liquid  composition,  if 

only  the  vapor  composition  has  been  measured   36-38  •   Finally,  if 

all  four  properties  are  known,  then  the  thermodynamic  consistency 
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of  the  results  may  be  tested  with  the  Duhem-Margules  equation  [39] 

The  most  widely  used  test  for  consistency  is  that  developed 
by  Coulson  &  Herington  [^o]  and  by  Redlich  &  Kister  Ul]  .   If  the 
volatility  ratio  is  defined  by  (X =  x^/xj  /  x^/x^  where  B  is  the 
more  volatile  component,  then  the  difference  in  excess  chemical 
potentials  of  two  components  can  be  given  as 

#B  -•**!  =  HT  ^T^Jj,)1   =  (£>CE/dxl)T  =  RT  ln(tXPA/PB)+(P-PA) 

(V1  "  V  "  (P  "  PB><VB  -  BBB»  +  2P"5B)AB(xf-  xf )  (k?) 
By  integrating  this  expression  over  x1,  where  G5  is  zero  at  x=0 

D 

and  x=l,  we  obtain 

£»(«  PA/PB)dxB  =  1/RT_['[(P  -  PB,(V1  -  BBB)  -  (P  -  PA)(vA  -  Baa) 

-2P(6B)AB(xf-  xf  )]dxB  (48) 

If  the  results  are  consistent,  the  volatility  ratio  will  equal 

the  pure  component  pressure  ratio,  ie.cx  =Pn/P»,  at  the  comcosi- 

E 
tion  where  G   is  a  maximum  or  a  minimum  value  in  a  non-ideal  mix- 
ture in  equilibrium  with  a  perfect  gas.   Therefore,  by  plotting 
^PA/PB  as  a  function  of  xB  the  areas  above  and  below  the  axis 
will  cancel  if  the  results  are  consistent. 

Unfortunately,  it  has  been  demonstrated  [28  J  that  this  test 
fails  for  liquid  mixtures  in  the  vicinity  of  their  normal  boiling 
points,  because  it  assumes  that  the  vapor  behaves  as  a  perfect 

gas.   Rowlinson  [lj  modified  the  method  by  noting  that  ((5b) 

AB 
made  no  significant  contribution  and.  that  P  could  be  replaced  by 

Raoult's  law  pressure,  x*PA  +  x^Pg,  so  that 
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|ln(«  PA/PB>dXB  =  [(PA  -  PB)(VA  -  VA  -  Baa  .  bbb)]  /2RT     (49) 

an  explicit  relation  which  is  convenient  to  use. 

The  most  widely  noted  method  of  determining  GE  from  measured 
isothermal  values  of  total  pressure  and  liquid  composition  is  that 
of  Barker  [^l\.      Th*s  author  has  made  calculations  of  GE  for 
N2"°2  mixtures  at  intervals  of  10°K  from  80°K  to  120°K  using 
this  method.   The  procedure  is  outlined  for  interest  in  Appendix  B. 

Rowlinson  [lj  has  indicated  the  difficulty  in  defining  the 
excess  free  energy  at  vapor  pressures  greater  than  2  to  3  atm, 
because  of  compressibility  effects.   This  can  be  seen  in  Fig.  32 
where  the  values  of  GE  calculated  by  the  author  have  been  plotted 
with  those  calculated  by  other  investigators  at  lower  tempera- 
tures.  It  should  be  noted,  however,  that  the  slope  of  the  values 
between  100  and  120°K  (dashed  line)  is  identical  to  that  at  lower 
temperatures  (solid  line),  and  one  wonders  if  perhaps  a  method 
for  calculating  GE  including  the  third  virial  coefficient  might 
not  displace  these  values  into  line. 

At  any  rate,  studies  in  the  vicinity  of  the  two-phase  region 
at  pressures  higher  than  atmospheric  are  virtually  non-existant 
[lt^3,*l4 J.   Even  the  molecular  theorists  have  contented  themselves 
with  predictions  of  excess  functions  in  the  liquid  phase  at  neg- 
ligible pressures  or  when  venturing  into  the  field  have  confined 
themselves  to  predictions  of  relations  for  the  critical  constants. 
In  the  gas  phase  preoccupation  has  been  with  determining  an  equa- 
tion of  state  [^5j  and  almost  no  experimental  work  has  been  done 
to  determine  excess  properties  "of  gases.   Only  very  recently  have 
badly  needed  measurements  on  dense  gases  been  attempted  ["^6,^7"]. 
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It  had  been  generally  felt 


UUj  that  the  molecular  theory 


would  stand  up  well  in  the  estimation  of  excess  functions  of  mix- 
tures of  condensed  gases,  that  is,  liquid  mixtures  of  simple  mol- 
ecules at  temperatures  near  the  boiling  points  of  the  mixtures. 
In  the  past  decade,  however,  experimental  work  by  Beenakker  et  al 

2i3»M  and  Staveley  et  al  I  6 , 7 J  have  shown  serious  inconsisten- 
cies and  plaguing  weaknesses  in  the  molecular  methods.   They  have 
•found  that  values  of  excess  energy  calculated  by  the  molecular 
theory  of  Prigogine  et  al  I  ^3J  are  strongly  dependent  on  both 
the  set  of  potential  parameters  chosen  and  on  the  pure  component 
selected  as  a  "reference"  substance. 

The  most  probable  cause  of  failure  of  the  theory  is  consi- 
dered to  be  the  Lennard-Jones  potential  |_  6,7 ,^-8,^9 1  50 J  particular- 
ly at  the  close  molecular  distances  involved  in  liquids.   Still 
another  cause  may  lie  in  the  innaccuracy  of  the  interaction  par- 
ameter combining  rules.   At  any  rate,  the  spectacular  differences 
in  the  value  of  excess  functions  calculated  from  molecular  theory 
from  those  measured  experimentally  is  enough  to  preclude  at  this  time 
consideration  of  values  formulated  by -this  means. 
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Calorlmetry 

The  measurement  of  the  excess  .heat  is  of  such  importance  as 
to  deserve  a  separate  section  reviewing  devices  and  methods  pres- 
ently available  for  its  determination.   The  two  principal  instru- 
ments are  the  adiabatic  and  isothermal  calorimeters.   The  adia- 
batic  calorimeter  reflects  positive- excess  heat,  if  the  tempera-- 
ture  drops  on  mixing  two  pure  components  in  an  adiabatic  contain- 
er.  If  on  mixing  or  afterwards  heat  is  added  to  cancel  out  the 
termperature  drop,  the  amount  of  excess  heat  can  be  measured. 
Heat  added  during  mixing  reduces  corrections  necessary  for  heat 
transfer  to  the  environment.   A  similar  determination  may  be  made 
for  negative  excess  heat  with  its  indicative  temperature  rise. 

The  adiabatic  calorimeter  has  been  bothered  by  a  variety  of 
problems  including  imperfect  heat  isolation,  incomplete  mixing, 
inaccurate  measurement  devices,  and  spaces  in  which  vapor  can 
form.   The  last  problem  is  complicated  by  the  need  to  allow  for 
volume  changes  on  mixing. 

Isothermal  calorimeters  have  the  advantage  of  not  requiring 
corrections  for  heat  losses  or  gains  and  are  able  to  measure  a 
greater  temperature  range  accurately. 

Many  such  devices  are  mentioned  in  the  literature  and  some 

of  the  more  important  are  listed  here.   Some  early  calorimetry 

•o 
at  temperatures  of  about  100  K  was  done  by  Rowlinson  &  Malcolm 

111.   A  versatile  instrument  which  can  be  used  either  isothermally 

or  adiabatically  has  been  developed  by  Nicholson  2 I.   A  glass, 

highly  accurate,  apparatus  is  described  by  Larkin  &  McGlashan   3  , 

which  has  eliminated  vapor  spaces,  provided  complete  mixing  and 
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allowed  for  expansion  and  contraction  through  a  compensating 
piston.   Likewise,  Mrazek  &  Van  Ness   ':   in  their  equipment  con- 
trol the  volume  with  close  clearance  pistons.   Another  novel  con- 
trol on  volume  change  is  the  metal  aneroid  diaphragms  of  the  de- 
vice of  Brown  &  Fock   5  •   Work  has  been  done  with  Bunsen  iso- 
thermal calorimeters  by  some  investigators  [6,?J,   Two  recent 
calorimeters  were  specifically  designed  for  measurements  of  con- 
densed gas  mixtures   8,9  . 
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Latent  Keat 

The  latent  heats  of  vaporization  for  mixtures  are  covered  by  a 
multitude  of  definitions   1,2  ,   The  most  obvious  are  the  heats 
required  for  vaporization  or  condensation  at  constant  pressure, 
constant  temperature  or  constant  composition. 

One  method  of  calculating  the  latent  heat  is  by  employing 
variations  of  the  Clausius-Clapeyron  relation   1,3-8  .   One  of 
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these  variations  expresses  the  differential  heat  of  vaporization 

AvvidP  -   (AHvi/T>dT  -  (*'  -  xlH  d?Gl/^x1?)VtTdx1  =  o  (50) 

for  a  snail  amount  of  vapor  phase  out  of  or  into  a  large  amount 
of  liquid.   Similarly,  the  differential  heat  of  condensation  is 

AvlvdP  -  (AHlv/T)dT  -  (x1  -  xv)(d2Gv/5x12)pjTdxv  =  0     (5D 

for  a  small  amount  of  liquid  phase  out  of  or  into  a  large  amount 
of  vapor.   The  volume  change  is  given  by 

Avvl  =  vv  -  v1  -   (xv  -  x1)(  dvVdx1)^  (52) 

and   the   heats  are 

AH        =   Hv    -   H1    -    (xv    -   x1)(dH1/dx1)p  (53) 

and  represent  the  heat  required  at  constant  temperature  and  pres- 
sure to  obtain  one  mole  of  vapor  phase  from  a  large  quantity  of 
liquid  or  to  condense  one  mole  of  vapor  into  a  large  amount  of 
liquid.   Similar  expressions  apply  for  /\  V   and  AK   .   Even  with 
this  seemingly  straightforward  approach  there  is  no  simple  method 
of  solution  except  with  assumptions  which  ignore  the  effects  of 
non  ideality  (see  Graphical  Methods), 

Corresponding  states  theory  yields  an  imposing  array  of  esti- 
mation techniques  for  pure  component  latent  heats  I  7  I  •   Calcula- 
tions were  made  by  one  of  these  methods  [9  ,  alleged  to  be  among 
the  most  accurate  and  simple  to  apply,  by  adapting  it  to  mixtures 
by  use  of  the  pseudocritical  constants.   This  empirically  derived 
relation  can  be  exoressed 
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<A;Wb'=  ['V5  log  Pcm  -  2.17)]  /(0.93  -  Tbr)  (5*0 

where  T^_  is  the  normal  boiling  point  and  A  H  ,  is  the  isotherrral 
b  «-j  vl 

differential  heat  of  vaporization.   The  results,  calculated  at 
one  atm,  were  less  than  satisfactory  as  can  be  seen  in  Fig.  3L 
In  fact,  because  of  the  large  latent  heats  as  compared  to  the  ex- 
cess heats  involved,  it  is  a  more  accurate  estimate  simply  to  use 
the  mole  fraction  average  of  the  pure  component  latent  heats, 
particularly  at  low  pressures. 

One   experimental  determination  of  the  latent  heats  of  N2-O2 
mixtures  exist   lOj ,  a  notable  number  when  so  few  for  any  binary 
mixture  can  be  found.   Dana's  measurements  at  atmospheric  pressure 
are  plotted  in  Fig.  31.   They  represent  the  integral  heat  of  vap- 
orization at  constant  pressure,  that  is,  the  heat  required  to 
change  one  mole  of  liquid  to  the  vapor  state  under  conditions  of 
constant  pressure  and  composition. 

No  matter  how  the  heat  of  vaporization  is  determined,  how- 
ever, it  is  only  as  useful  as  the  Isotherm  or  phase  boundary  from 
which  it  will  be  based. 
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Graphical  "ethods 

The  earliest  graphical  techniques  for  determining  the  thermo- 
dynamic properties  of  N?-0   mixtures  were  devised  at  the  Kamer- 
lingh  Onnes  Laboratory  at  the  University  of  Leiden  I  1 v  2 J .   Keesom 
&  Tuyn  developed  a  general  construction  for  HtS,x  surfaces  of  the 
mixture  in  the  two-phase  region  at  atmospheric  pressure.   The 
property  of  the  mixture  that  is  the  key  to  the  geometric  method 
is  the  transverse  fold  in  the  H, S ,x-surface  representing  the  wet 
vapor  region.   A  tangent  plane  across  the  fold  contacts  the  points 
of  equilibrium  or  coexistant  phases.   These  points  were  called 
connodes;  the  isotherm  joining  them,  the  nodal  line;  and  the 
locus  of  the"se  points,  the  connodal  curves  (vapor  and  liquid 
curves ) . 

Basically,  their  procedure  was  to  locate  known  points,  such 
as  the  boiling  point  of  oxygen  on  an  H-S  diagram,  and  with  the 
aid  of  pure  component  heats  of  vaporization  at  l'atm,  the  defi- 
nition of  enthalpy  and  entropy  in  terms* of  specific  heat,  ideal 
mixture  rules,  and  Dodge  and  Dunbar  equilibrium  data   3  »  plot 
the  vapor  line  which  is  then  easily  projected  into  the  S-x  plane. 

Using  the  previously  mentioned  principle  of  tangent  planes. 
the  liquid  line  was  determined  in  the  H-S  and  S-x  planes  by 
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graphically  constructing  the  line  representing  the  plane  tangent 
to  the  connodal  points,  assisted  by  Dana's  heats  of  vaporization 


^ 


The  projection  in  the  H-x  plane  was 


at  constant  composition 
then  easily  established. 

The  procedure  then  recommends  extrapolation  in  the  H-S  plane 
of  the  lines  of  constant  composition  to  higher  pressure  using 
Second  Law  relationships  and  the  vlrial  equation  of  state,  ie, 


dH  =  V  -  T(dv/$T)     dP    and 

r  ,  X 
dS  =  -(c)V/c)T)PfXdP 


which  become  with  the  application  of  the  equation  of  state 


H  =  (Pg  -  P2)(B  -  TdB/dT)    and 


S  =  R  In  P1/P2   -  (p2  -  pi)  dB/dT 


M55) 


f  (56) 


The  ideal  mixing  rules  locate  the  isocomposition  lines  between 
those  of  the  pure  components,  and  the  equilibrium  data. fixes  the 
end  points  at  the  vapor  line.   The  form  of  the  liquid  line  is 
then  taken  to  be  the  same  as  that  for  one  atm. 

Lerberghe   5   proposed  a  variation  on  this  method  of  his 
colleagues  which  requires  only  the  use  of  the  T-x  equilibrium 
diagrams.   From  thermodynamics  and  analystical  geometry,  expres- 
sions were  developed  for  the  enthalpy  and  entropy  difference  be- 
tween the  connodal  points. 


Hv  -  H1  =  (N  +  NT)(xV  -  x1) 
Sv  -  S1  =  (L  +  M)(xV  -  x1) 


and 


>  (57) 
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whe 


v         v  r,  v .  ,  vi 

re      L   =   SA   -   S3  +   R    In  1(1    -   x    )/x 

M  =    -RT[l/(l-xV)    +    l/(xV)]dxV/dT 


>   (58) 


where  the  A  and  B  are  the  pure  component  properties  at  the  pres- 
sure and  temperature  of  the  mixture,  extrapolated  as  necessary 
into  the  wet  vapor  region. 

Thus  by  fixing  the  vapor  line  as  was  done  by  Keesom  and  Tuyn, 
Lerberghe's  method  locates  the  liquid  line  with  the  help  of  the 

equilibrium  data.   Note  that  the  requirement  exists  to  perform 

v 
the  difficult  graphical-  tangent  determination  for  dx  /dT  or  to 

derive  an  equation  for  the  vapor  line  on  the  T-x  diagram  by  curve 

fitting  and  then  taking  its  derivative  analytically,   Lerberghefs 

graphical  measurements  are  reflected  in  Fig,  31  for  the  liquid 

line  at  one  a  tin. 

Several  years  later  Bosnjakovic   6  nicely  summarized  vari- 
ous graphical  procedures.   His  method  suggests  that  in  most  cases 
the  excess  functions  in  the  vapor  phase  can  be  neglected,  at 
least  in  moderate  pressures,  and  the  isotherms  in  the  H-x  and  S-x 
diagrams  may  be  drawn  according  to  the  ideal  mixture  rules.   Then, 
if  the  temperature-vapor  composition  relationship  is  known  from 
equilibrium  data,  the  vapor  line  is  also  determined  in  both 
planes  0 

In  the  liquid  region,  if  the  excess  properties  are  known, 
then  the  isotherms  may  be  drawn  and  the  liquid  line  again  deter- 
mined as  above  for  the  vapor  line.   If  the  excess  heat  is  not 
known,  some  estimation  of  the  latent  heat  of  vaporization  must 


be  made.   The  Clausius-Clapeyron  relation  (Eq, -51)  makes  a  good 
starting  point.   From  the  definition  of  Gibbs  free  energy  the 
second  derivative  is  given  by 


2 


2 


(  6''G/0x%fT=  (  ()V5/)ptT-  T(  ^'S/5x2)P(T 


Then  for  an  ideal  vapor  mixture 


2 


(59) 


(dV/^xV2)    .  RT/xV(l-xV) 


(60) 


Measurement  of  the  slope  of  the  T-x  diagram  graphically  or  analy- 
tically yields  (  d  x  /dT)p«  Enough  information  is  then-  available 
to  solve  Eq.  (51)  for  /\h,   at  constant  pressure  and  temperature. 

The  vapor  isotherm  is  then  extended  to  intersect  x  and  /\  H 

lv 

is  plotted  from  that  point  (Fig.  5).   Repeating  this  process 
with  isotherms  of  the  wet  vapor  region  yields  the  liquid  line 
(Fig.  31).   Obviously  the  method  is  sensitive  to  the  behavior  of 
the  two  differential  terms  of  Eq.  (51), 


AHlv 


'J- 2 


In  the  case  of  the  S-x  diagram,  .after  the  isotherms  have  been 
drawn  according  to  the  ideal  mixing  rules  and  the  vapor  line  de- 
termined from  equilibrium  data,  the  liquid  line  may  be  determined 
from  the  relation 

Ah1v/t  =  s1  -  sv  -  (x1  -  xv)(5sv/5xv)pj  (61) 

The  value  of  /\  H,  /T  can  be  obtained  from  the  H-x  diagram  and 
plotted  from  the  intersection  of  the  tangent  to  the  vapor  line  at 
x  with  the  equilibrium  liquid  composition,  x   (Fig.  6). 

A  secondary  method  of  constructing  the  S-x  diagram  when  the 
H-x  diagram  is  available  utilizes  an  integration  procedure,  which, 
needless  to  say,  is  less  prone  to  graphical  error  than  the  dif- 
ferentiation technique.   This  method  uses  the  relation  that,  for 
a  reversible  vaporization  process  at  constant  pressure  and  con- 
stant liquid  composition, 


-  SV  -  J,1 


(SV  -  Sx)_l  =   /  1/T  (dH)  1  (62) 

i         x 
T1 

The  right  hand  side  of  the  equation  may  be  easily  integrated 
graphically  using  any  of  the  several  numerical  integration  rules 
and  the  difference  measured  from  the  known  value  of  S  , 
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The  Thermodynamic  Properties 

The  thermodynamic  properties  of  binary  mixtures  of  nitrogen 
and  oxygen  considered  to  be  of  the  most  value  for  use  in  engineer- 
ing calculations  are  those  of  enthalpy  and  entropy.   A  useful 
format  for  the  presentation  of  these  properties  is  the  graphical 
display,  in  particular  the  enthalpy-composition  and  entropy- 
composition  diagrams  at  constant  pressure.   Accordingly  Figs, 
7-30  have  been  prepared  for  mixtures  at  pressures  of  1 ,  5#  10, 
and  20  atm.   Enthalpy  and  entropy  have  been  represented  as  func- 
tions of  composition  in  the  liquid  line  region,  the  vapor  line 
region,  and  the  wet  vapor  region. 

Values  for  the  pure  component  oxygen  were  taken  from  Stew- 
art's machine  tabulation   1   based  on  a  complex  empirical  expan- 
sion.  Those  for  nitrogen  came  from  the  work  of  Strobridge  I2J 
whose  machine  tabulation  is  likewise  based  on  an  empirical  equa- 
tion of  state. 

The  molecular  weight  of  nitrogen  was  taken  as  28.016 
grams/mole   2  and  one  calorie  as  il-,18^  abs ,  joules,  by  which 
conversion  factors  the  tabulated  values  of  Strobridge  were  con- 
verted to  abs.  joules/mole. 

For  the  sake  of  graphical  convenience  the  properties  of 

o 

nitrogen  and  oxygen  were  equalized  at  91  K  at  1  atm,  creating  an 

oxygen  enthalpy  scale  factor  of  ^2^5.3  j/mole  and  a  nitrogen  en- 
tropy scale  factor  of  81.12  j/mole.   The  values  of  entropy  are 


Ift 


relative  only,  since  the  nitrogen  scale  is  not  based  on  zero  en- 
tropy at  absolute .zero  as  is  that  of  oxygen.   Property  values  of 
the  pure  components  within  the  two  phase  region  have  been  linear- 
ly extrapolated. 

Values  of  the  entropy  of  mixing  (Gibbs  paradox)  are  derived 
from  the  relation 


,ra 


/.x.  lnx. 
*   i     i 


S   =  R  Z_iX,  lnx,  (63) 

_2 

where  R  is  the  universal  gas  constant  equal  to  8.2057^-  x  10' 

L-atm/mole-°K  [l] . 

Isotherms  in  the  superheated  vapor  region  were  derived  as 
mole  fraction  averages  of  the  values  for  the  pure  components,  ie, 
from  the  ideal  mixture  rules.   The  excess  heat  and  entropy  func- 
tions for  the  liquid  isotherms  were  derived  by  the  author   from 
the  data  of  several  investigators  (see  Discussion  of  Results). 
The  vapor  and  liquid  lines  were  plotted  using  the  isotherms  and 
the  temperature-composition  data  of  Dodge  and  Dunbar  |3|« 

References  for  Thermodynamic  Properties 

la   Stewart,  Ph.D.  Thesis,  Univ.  of  Iowa,  1966. 

2.  Strobridge,  MBS  Tech  Note  129,  1962. 

3.  Dodge  &  Dunbar,  J.  Am.  Chem.  Soc.,  ^9,  591,  1927. 


45 


r\4 


66-- 


58-- 


50 


x  42 

UJ 

o 


>- 


-J  26- 


< 
x 

-z. 

UJ 


18  - 


10- 


2- 


-p                                            -  

/ 

• 

90°K 

/88K 

/86 

K     /84K 

/82K  /80K 
/          /     1 

•  / 

/  78°K 

p/"  i/"ml 

i 

Wt  i   /v/-\ruK/K 
N2-02/  1  ATM/ 

EGION 

/ 

, 

t 

/    ' 



i 

r^< 

^1 

i 

10      20     30     40      50     60,    70     80 
COMPOSITION    %M0LE   N2 

Figure    7 


90    100 


46 


1 

i 

[ 

1 

i 

TV 

1^- 

00 

o" 

< 

DO 
30 

c 

>CD 
00 

00 

/  * 

oc 

i 

• 

Ul 

z 
Ij 

O 

Oh 

l 

>  z:  - 

i 

c 

i 

1          \               1           II 

o 

00                  N-                   CD 

If)            55- 

ro 

cd 

CO                CD                 CD 

CD               CD 

CD 

pi  x  310 w/r 

AdlVHlNB 

o 

o 


o 
o 


o 

00 


o 


o 

CD 


O 
C\J 


O 


C4 


o  uj 

If)  -J 

o 


o 


o 
o  — 

CO 
O 
Q_ 


O 


Figure  8 


47 


. 

\ 

_ 

v 

1 

1 

1 

1 

\ 

s. 

LU 

—I 

\ 

\ 

v^       • 

Q 

_J  2 

\ 

v^ 

\ 

x^ 

2C 

°00 

SsN- 

°o 

N^ 

■ 

00 

. 

-  co 

-  co- 

• 

I 

1 

• 

i 

°o 
o 

1 

1 

1 

CO 


o 


m 


ro 


CM 


01  X310W/P      Ad1VHlN3 

Figure     9 


O 

o 


o 


O 
CO 


n-z: 

LU 

_J 
oo 

.o 
o^ 


o 

in 


o 


o 


(j) 


a. 


o 


O 


O 

CM 


O 


O 


48 


UJ 


LU 


170 


160- 


150- 


140- 


130- 


o 
a: 


0 


100 


90^ 
0 


I 

i     i 

^= 

' 

/j 

i^-^^ 

— ~r 

• 

"/ 

•      1 

90°K 

/88°K       I 

8S°K 

/84K    / 

82°K 

/80K     1 



/      78°K 

_  I 

. 

1  WE' 

r  w 

^POR 

/REC 

ilON 

NfOgh  ATM    / 

1 

i 

f  • 

* 

yy 

1 

• 

10      20      30     40      50    .60     70      80      90     100 
COMPOSITION     %M0LE    Na 

Figure   10     • 


49 


O 

o 

— 

/ 

"'"'  A>fjrsts'i^ 

s%/yy^ 

y^ 

I 

CO 

o 

/Y /\/v  v 

G) 

1* 

00 

O 

-1 

'    A 

CO 

o 

■]—/— 

CO 

oo 

z 

UJ 

^O 
O 

O 

< 



« 

CO 

0          5 
S1TI0N 

v  _ 

0          4 

COMPC 

\ 

.    loo 

UJ 

z 
-J 

oc 

ro 

O 

\ 

5  r5 

c\J 

o 

°o 

CL    ,    •— 
<    U< 
>Z  — 

1 

1    1 

I  N. 

O 

s 
r^ 

c 

c 

3 

310 

w/r 

Pigi 

:> 
A 

are 

dOcJ 

ii 

IN] 

3 

c 

J 
D 

50 


-g1 

CO 

1 

1 

1 

1 

y^. 

^n 

1 

1 

i 

"If 

LU 

• 

• 

• 

-J 

_co_ 

Q     »  ^ 

<* 

-J    z    - 

X 

<- 

CO 

CO 

l\ 

• 

-  en  - 

°co 

CO 

. 

eft 

1 

1 

1 

1^^ 

^r^- 

^^ 

1         1         1         1         1         1 

■'■  N 

-^ 

2s». 

^^J 

^ 

"^ 

^i 

o 
o 


o 

CD 

O 

CO 

<nj 

O 

z: 

LU 

-J 

o 

o 

ID 

2> 

^o 

o"" 

o 

Z 

o 

1— 

o 

en 

^r 

o 

CL 

s 

o 

o 

u 

ro 

o 

CM 

o 


o 

CO 

o           sr 

c\I               O 

CO 

o 

CD 

CD                 CD 

CD                CD 

co 

>L310W/r 

Ad0UlN3 

Figure   12 


51 


Figure   13 


70 


62 


54 


C\j 


X 

LU 

o 


46  — 


*38 


< 

x 


lu22 


14- 


6$ 


, 

1    1    1    1 

'               1 

1         1 

1                                    |                                      . 

i 

— 

T"6^"" 

I08K 

llOGK    /I04K 

/|020k/|00k/98K    |96K 

1 



—  1 



—  /  — 

/  WE 

T /VAPOR  REGION 

/N2 

UzpAl 

Iwl 

i 

1 

1 1 

I 

10 


20     30     40     50     60     70     80 
COMPOSITION     %M0LE    N2 


90    100 


52 


1 

1 

1 

I 

/ 

1 

1 

o 
CO 

i 

D 

o 
o 

/ 

o 

r 
C 

c 

1   / 

{ 



o 

o 

— — / 



6 
< 

< 

o 

co 
o 

i 

30  — 

_j 

t 

1 

>  Z 

H 

i  < 
IT) 

l 

O 

O 

co 

h- 

CO 

if) 

r^- 

CO 

CO 

CO 

CO 

CO 

01  X  310W/P        AdlVHl 


N3 


o 
o 


o 


o 

CO 

Z 

O 

l^- 

LU 

-J 

o 

2 

o 

o 

CO 

o 

o 

Z 

m 

o 

h- 

o 

10 

sr 

o 

Gl 

^ 

o 

o 

u 

ro 


O 


Figure   14 


53 


s 

\ 

i 

1 

1 

1 

1 

1 

0 

QJ 

z 

Q 

a  i 

-J   2 

\ 

\ 

<   2 

:  in 



°CO 

r 

> 

• 

V 

o* 

• 

. 

'  c\ 

c 

J 

) 



* 

C 

<• 

\ 

\ 

—  * 

c 

*» 

\ 

5 

> 

V^v^ 



V 

c 

s 

3 

iii             i             i 

i   i   i   i   i   i   i 

1                 '                 I                 !                i 

1        1        1        1        1 

1 

\^ 

o 
o 


o 


o 

co 


o  z: 


O 
If) 


o 


O  </) 

KT  O 

Ol 


ro 


O 


O 


J o 


o 


lO 


rO 


oj 


O 


B- 


01  X  310W/P      Adl VH1N3 


Figure   15 


54 


160 


150 


o 
LU 


-i  140  — 


o 

21 


30 


Q. 

o  120 

z 

UJ 


110- 


100^ 


i 

1    1    1 

1       ' 

7^ 

i 



D8K 

/I06H 

;  /io< 

4K/IC 

0          i 

)2K  i 

I00°K 

/   ° 

/98K 

;    k 

0 

SK— 

i 

r^ 

_^Vs^ 

■    "    i 

.  t     ir  t      \/Ar 

1         1 

WET  VAF 
N2-0a. 

^OR 

5  ATM 

10 


20 

30 

40     50 

60 

70 

80 

COMPOSITION 

%M0LE 

N2 

Figure 

16 

90    100 


55 


/ 

X^i 

~~1 

/ 

/// 

00 

/  / 

/ 

o 

°c\j 
o 

o 

o 



— \ — 

o 

LU 

-O    -    S 

— — 

O 

Cl 

■z  u 

c 
"> 

CD 

sr 

<\J 

O 

co 

CD 

sa- 

CD 

CD 

CD 

CD 

m 

in 

in 

>!03"10W/P       AdOUlN3 

Figure   17 


94 
m 


o 
o 


o 


o 

1^- 


o 

CD 


LU 

o 


o 
m 


o  en 
a. 


o 


O 


O 


O 


o 


56 


1 

1     1 

/\ 

1 

— 1 

i 

i 

LU 

-J 
Q 

|98K 

3     <M  2 

CfOh 

-I    Z    ID 

o 

o 

\ 

o 

I 

V 

o 

\V 

I08K 

w  \  N{\  7\  fSv 

^b^H^ 

l^s^ 

1    1    1    1    1    1    1    l> 

1  \ 

l\^ 

i\ 

c\J 


CO 

o 


ID 
O 


O 


CM 
O 


O 
O 


>!   3"10W/r  Ad0dlN3 

o  ' 

Figure   18 


O 
O 


o 


o 

CO 


—  ~    ot 


LU 


O 

CO 


O 

m 


a. 


O 


O 


O 


57 


70 


62 


CO 

»0 


54 


x 

w    46 
o 


38 


>- 

Q. 

-J 

< 
X 
I- 

z 
ui 


30 


22 


14^ 


1 

1                 1        1 

— 

II8°K 

II  I6K/1I- 

OK/1 

2°k/| 

I0°K 

|i08°H 

:    I06K 

1     I04°K 

_J 

^^~-~ 

■p^^ 

h— WET    V 
N2-02 

APOf 

3    RE 

Gior 
oa™ 

j — 

4 

0 


20     30     40     50 
COMPOSITION 


60     70     80 
%MOLE    Hz 


90    100 


Figure   19 


50 


'| 

II 

r~r 

1 

i 

il 

II  1 

t   1           1 

/ 

1 04  K 

1 

— -.-=: 

^1 

°00 

°1 

:] 

V 

^1 

°o 

V" 

1/ 

y 



> 
o" 

C 

LU 
i    _ 

OC 

, 

O 

CL    < 

2 

3   < 

1     ^ 



(\J 

>l     / 

>  z   - 

• 

i 

ro 


c\i 


O 


o 
cd 


CD 
CD 


CD 


CD 
CD 


O 
O 


O 
(J) 


O 

CO 


O 


M 


UJ 


O 
CD 


O  2; 

L00 


O  CO 
CL 


O 
cvj 


_o 


,01  X310W/r       AdlVHlN3 


Figure   20 


59 


ro 

CM 


CM 

cvj 


O 


—  O 

CM  CM 

,01    X  310H/f  Ad1VHlN3 


oo 


Figure   21 


60 


20     30     40 
C0MP0SIT 


50 
ON 


60      70     80 
%MOLE    Hz 


90    100 


Figure  22 


61 


1 

1  j^y 

/aA// 

x*- — 

O 

i — 

,.,,., 

i —  ■ 

/ 

/ 

'    / 

/ 

-I06K 

1 

/ / 

CO 
O 

/ 

—  -  - 

/ 

' 

o 

i 

u. 

z 

J 

1 



o 
CO 

c 

\       CM    -2. 

:    i  •<- 

:  ^2 
1 

< 
> 

^ 

\* 

N, 

"°co" 

7 

1 

o 

oo 

CO 

^r 

r\] 

O 

00 

CO 

m 

m 

m 

lO 

in 

<sT 

>i0310W/r  Ad0tiiN3 

Figure  23 


O 

o 


o 


o 

00 


2  - 


o 

CO 


o 
m 


LU 

o 

^o 
o 


O 
I- 


O    CO 

sr  o 
a. 

^  ° 

8" 


O 


O 


62 


1 

o 

i  j 

^1 

• 

r  i 

1 

V 

/[/ 

o 

« 

UJ 

z 

00 

.  _J 

o 

Q 

ID 

_J 

2: 
0  ^ 

"1 

1  * 

2?< 

t 

3 

°§ 

— 

« 

— 

°S 

•V 

— 

•£ 

[""v     [S. 

nK 

1 

^^ 

xr^ix.r^ 

^o\ 

\ 

1  XZr^TNw  >v>^rvr — 

r^^ 

<D 


C\J 


O 


00 
O 


o 


o 


>l  310W/r  Ad0UlN3 

Figure   24 


O 
O 


o 
0 


o 

00 


UJ 

o  -» 
iO  o 


o 


o 


sr  o 

CL 


o 

o  u 

ro 


O 

cm 


O 


O 


63 


X 

LU 

o 


70  — 


62 


c\,      54 

'o 


46 


38 


<• 

x  3( 

h- 
z 

LU 


22^ 


I       ! 

1 

1 

1 

i 

i 

i 

• i 

1 

- 

ir™ 

:=r 7* ~ 



• 
30K 

32°k|i: 

28K 

fl26°h 

;/|24k/|22°k/|2( 

D°K  1 1 

1  1 

0 

8K 

0 

6K 





WE" 

r  va 

POR 

REGIC 

-Nft-02    - 

20  ATM 

)N 



14  — 


0       10 


20     30     40     50     60    70     80     90    100 
COMPOSITION      %M0LE   N2 


Figure   25 


64 


CO 


O 


o 

Is-  o 

0!  X  3iow/r 


CD 


cd 


CD 
CD 


AdlVHlNB 


Figure  2  6 


65 


\ 

1 

i 

I 

j 

"1 

i 

LU 

-J 
Q 

1 

i 

\ 

■i 

°i 

OJ 

51 
*>    1- 

^s_ 

V 

— 1 

\ 

-J 

z:  oj 

\ 

oj 

|\ 

\. 

N^- 

CO 
OJ 

\ 

I^S^ 

ro 

\ 

"^^ 



OJ 

ro 

\ 

V 

111                 i                         ii                                    ""V* 

i    i    I    I    1    i    i    i  ■  I    I    I    I    i    I  > 

\     ,  . 

ro 


ro 


o 

G> 

oo 

r^- 

<s 

ro 

OJ 

OJ 

OJ 

OJ 

,OIX3H0H/r       AdlVHlN3 

Figure   27 


O 
O 


o 
o 


o 

CO 


O    2? 


LU 

_J 
o° 


O 
IT) 


O 


O  CD 

sr  o 

Q. 

o 


o 

ro 


O 
oj 


o 


J O 


66 


0       10      20     30     40     50     60     70      80     90    100 
COMPOSITION     %MOLE    N2 


Figure   28 


67 


O 

j*  i 

°CO 

1 

1 

O 

o 
o 

// 

y     \ 

A 

°o 

<\J 

/      /       / 

00 

1 

^ 
M 

o£ 

I  r 

LU 

o 

9^ 

,     __ 

\    \ 

cm 

^     - 

'to 

CM 

&    o 

~~ 

t 

DL 

• 

•J 

0 

- 

ro 

1 

UJ 

ro  w 
O 

CVJ 

O 



( 

« 

a 

c 

a 
< 

>  o  H 
.   1    < 

•     d  o 
:  z:  oj 

,rO 

• 

> 

r^ 

I 

1       1 

■  O 

{ 

I 

D 

u 

r 

C 

L 

n 

< 

D 
n 

c 

D 

-4- 

* 

o 

S 

S 

t 

T 

>{o310W/r  AdOdlNB 

Figure  29 


68 


•       1        1                                  / 

1 

/  — 

1 

o 
CO 

• 

<M 

r\  1  - 

LU 

Q 

o  ; 
i 

c<   < 

2     t 

£ 

OJ 

ID 
-J 

D 
M 

CM 

cm 

v: 

°Q0 
CVJ 

•g- 

ro 

\\ 

\. 





*cm 
ro 

CjvInJv 

(MSlXN 

\ 

\ 

1    1    1    1    1    1    1 

^^. 

^ 

\w 

\ 

OsJ 


CM 
CVJ 


o 

CM 


CO 


>{   310W/r  AdOUlN3 

0  ' 

Figure  30 


CM 


O 

O 


o 


o 

00 


o 


o 

CO 


-J 

o 


o 

O  a) 

a. 


ro 


O 


O 

CM 


O 


69 


DISCUSSION  OF  RESULTS 
The  choice  of  the  method  of  determination  and  the  form  of 
presentation  of  the  results  was  arrived  at  after  consideration  of 
all  the  previously  discussed  techniques.   No  attempt  was  made  to 
formulate  a  virial  expansion  or  other  empirical  equation  of 
state,  because  of  a  general  lack  of  the  PVT  data  required  for 
such  a  formulation.   A  corresponding  states  theory  estimation 
technique,  using  pseudocritical  constants  and  the  enthalpy  devi- 
ation tables  of  Lydersen.,  was  conducted,  with  the  results  already 
mentioned.   The  negative  excess  heats  calculated  by  this  method 
are  thermodynamically  inconsistent,  and  this  procedure  was  re- 
jected.  Using  pseudocritical  constants,  a  modified  relation  for 
the  latent  heat  was  calculated,  which  when  plotted  produced  the 
least  agreeable  of  all  liquid  lines  determined  (Fig.  31).   The 
purely  graphical  approach  of  Lerberghe  produced  better  results, 
except  at  high  molar  concentrations,  which  plotted  extremely 
erratically..  Both  a  mole  fraction  average  of  the  heats  of  vap- 
orization of  the  pure  components  and  a  mole  fract.ion  average  of 
the  pure  component  boiling  points  produced  better  results  in 
that  order,  than  did  the  complex  estimation  techniques  mentioned 
above.   A  plot  of  the  results  obtained  by  these  methods  at  1  atm 

together  with  the  method  finally  chosen  is  given  in  Fig.  31. 

> 

In  the  end  it  was  decided  that  the  calculation  of  excess 
functions  provided  the  best  determination  of  the  liquid  line,  and 
therefore  of  the  isotherms  in  the  liquid  region.   The  liquid 
line  determined  by  this  method  shows  very  good  correlation  with 
the  experimental  evidence  of  Dana,  as  can  be  seen  in  Fig.  31. 
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Figure  31.  Liquid  line  for  Ng-C^  mixtures  at  atmospheric  pressure. 


A  presentation  of  the  excess  Gibbs  free  energy  as  calcula- 
ted by  previous  investigators  is  shown  in  Fig.  J2,      The  excess 
Gibbs  free  energy  is  related  to  the  other  excess  functions  by 


F      F1       E      F        "\  E 

JT  =  G^  +  TS^  =  G"J  -  T(  OG  /£)T) 


(^6) 


.E 


It  is  clear  from  Fig.  32  that  the  slope  of  G  with  respect  to 
termperature  is  very  linear  in  the  low  pressure  region  and  has 
been  measured  as  0.08  cal/mole  °K.  With  this  slope  the  excess 
heat  of  the  liquid  was  calculated  to  be  16.0  cal/mole  at  equimolar 
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composition  and  low  pressure.   Only  one  experimental  calculation 
of  excess  heat  for  ^2-°2  ^-iclu^  mixtures  has  ever  been  made   1  . 
These  measurements  by  Knobler  et  al  of  the  excess  heat  at  77  K 
are  plotted  as  Fig,  33.   These  results  can  be  expressed  by  the 
fitted  polynomial  approximation  I 2J , 

HE  =  xN2  xQ2  [if 2  +  10(2x02  -  1)  +  2.0  (2xQ   -  l)2] cal/mole      (6*0 


which  yields  an  excess  heat  of  10.5  cal/mole  at  x=0.5.   A  concur- 

E 
rent  calculation  of  G  was  made  and  can  be  given  by 

GE  =  xN2  xQ2  (\o  +  M2x02  -  1)  +  ^(2x02  -  l)2]cal/mole.       (65) 

This  gives  an  excess  Gibbs  free  energy  of  10.0  cal/mole  at  x=0.5 

or  a  value  nearly  the  same  as  that  of  the  excess  enthalpy.   Since 

E  E   E 

the  temperature  derivative  of  G   is  negative,  then  H  :=>  G  ,  and 

from  the  weight  of  evidence  by  considerably  more  than  0,5  cal/mole. 

•It  was  decided  that  16.0  cal/mole  was  a  reasonable  value  for  the 

excess  heat  at  equimolar  composition  in  keeping  with  Fig.  32, 

and  the  good  fit  obtained  in  plotting  the  liquid  line;  Knobler' s 

data  was  then  scaled  up  to  this  result.   It  should  be  noted  that 

in  agreement  with  prediction   3   the  excess  heat  is  asymmetric 

With  respect  to  the  composition,  while  the  curve  of  excess  free 

energy  is  symmetrical.   Values  of  the  excess  entropy  were  then 

calculated  from  Eq.  (^6). 

When  the  calculated  values  of  excess  heat  were  applied  to 

isotherms  in  the  liquid  region  at  one  atmosphere,  the  resulting 

liquid  line  showed  excellent  agreement  with  the  experimental 

data  of  Dana".   It  was  decided,  therefore,  that  this  technique 
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Figure  J2,    Excess  Gibbs  free  energy  as  a  function  of  temperature 
for  equimolar  liquid  mixtures  of  N2-O2.  See  Appendix  C  for  data 
sources. 
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Figure  33.  Excess  heat  and  excess  Gibbs  free  energy  as  a  function 

E 
of  composition  from  the  data  of  Xnobler  et  al.  H  fitted  to  re- 
vised equimolar  avlue  of  this  thesis. 
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at  atmospheric  pressure  was  the  best  available. 

The  integration  method  of  Bosnjakovic  was  used  to  calculate 
(Sv  -  S  )xl.   The  liquid  line  thus  determined  for  the  S-x  dia- 
gram was  in  excellent  agreement  with  that  calculated  from  values 
of  the  excess  entropy. 

In  the  vapor  phase  it  had  been  decided  that  most  authors 
generally  agreed  that  there  was  negligible  heat  involved  with  - 
mixing  at  low  pressures  including  atmospheric,  so  that  the  best 
approximation  of  isotherms  was  a  simple  mole  fraction  average 
corresponding  to  the  mixing  rules  for  an  ideal  mixture;  this  in- 
cludes the  entropy  of  mixing  for  the  S-x  diagram  vapor  line. 

The  major  problem  then,  concerned  pressures  greater  than 
atmospheric,  that  strangely  ignored  and  abandoned  region  where 
investigators  have  feared  to  tread.   It  would  be  foolhardy  and 
certainly  not  justifiable  to  make  unwarranted  assumptions  about 
conditions  at  higher  temperatures  where  pressures  exceed- several 
atmospheres,  where  the  third  virial  coefficient,  the  compressi- 
bility of  the  liquid,  and  the  volume  change  on  mixing  are  all 
important  factors. 

For  the  best  available  presentation  it  was  decided  that  the 
excess  functions  of  heat  and  entropy  already  calculated  would  be 
used  to  determine  the  isotherms  of  the  liquid  and  mole  fraction 
averaging  the  vapor  Isotherms  at  pressures  and  temperatures  be- 
low the  critical.   This  decision  is  not  wholly  without  justifica- 
tion.  It  is  not  unreasonable  to  assume  that  the  excess  thermo- 
dynamic properties  of  the  liquid  are  relatively  independent  of 
pressure  [ifj ,  principally  because  of  the  relatively  small  volume 
changes  in  the  liquid  state.   Additionally,  Fig.  32  gives  some 
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small  hint  that  deviations  at  the  higher"  temperatures  associated 
with  the  higher  vapor  pressures  are  not  so  great  from  the  more 
accurately  calculated  linear  results  at  lower  temperatures. 

The  incompressibility  reasoning  does  not  apply  so  well  in 
the  dense  gas  region,  but  here  too  it  is  felt  that  the  ideality 
of  mixture  assumption  is  valid  enough  when  compared  to  other 
available  data.   The  values  of  several  investigators  for  air 

5»6,7,8  were  compared  at  the  pressures  in  question  at  several 
isotherms  near  the  dew  point  and  compared  to  the  79*mole  No  mix- 
ture of  nitrogen  and  oxygen  isotherms  calculated  by  mole  fraction 
averaging.   Although  the  deviations  varied  widely,  a  positive 
excess  heat  with  increasing  pressure  might  seem  to  be  indicated, 
the  effects  of  argon  in  the  mixture  and  the  particular  investi- 
gations notwithstanding.   The  deviations  are  not  large  enough  nor 
their  pattern  clear  enough  to  justify  guesswork  at  possible  val- 
ues of  excess  heat. 
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CONCLUSIONS 

The  intention  of  this  thesis  was  to  survey  the  literature 
with  regard  to  the  thermodynamic  properties  of  binary  mixtures  of 
nitrogen  and  oxygen  in  the  two-phase  region  and  to  analyze  the 
resulting  data,   Fromtthe  sizeable  collection  of  available  print- 
ed matter  which  concerns  itself  directly  and  indirectly  with  the 
subject,  it  would  seem  an  easy  matter  to  make  an  analytical  de- 
termination of  the  common  thermodynamic  properties  of  such  rela- 
tively simple  mixtures.   A  wide  variety  of  estimation  techniques 
do,  in  fact,  exist,  and  while  those  have  been  presented  here, 
most,  if ' not  all,  of  them  are  either  unsatisfactory  or  at  least 
only  partially  useful. 

Specifically,  the  virial  expansion  and  other  equations  of 
state  with  their  associated  derivatives  would  be  most  satisfac- 
tory if  experimental  PVT  data  existed  with  which  to  determine 
their  coefficients.   Kolecular  theory,  similarly ♦  has  great  po- 
tential, particularly  for  liquids  and  dense  gases,  but  suffers 
from  the  uncertainty  involved  in  its.  basic  parameters.   The  tech- 
niques involving  pseudocritical  constants  give  good  rough  esti- 
mates, but  do  not  provide  the  accuracy  required  for  even  engi- 
neering  calculations.   The  excess  functions  provide  a  rapid  and 
relatively  easy  method  for  determining  the  thermodynamic  proper- 
ties.  Presently  available  data  on   these  give  good  results  at 
pressures  near  atmospheric. 

All  of  the  aforementioned  methods  have  one  thing  in  common. 
None  of  them  have  been  developed  for  use  in  the  high  pressure 
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area  of  the  two -phase  region.   Where  no  one  has  gone  before,  how- 
ever,  any  determination,  no  matter  how  simplified,  idealized,  or 
approximated  It  may  be,  is  better  than  none  at  all.   Therefore  it 
was  decided  that  by  taking  the  excess  functions  independent  of 
pressure  in  the  liquid  region  and  using  ideal  mixture  rules  in  the 
vapor  phase,  satisfactory  thermodynamic  diagrams  could  be  con- 
structed. 

RECOMMENDATIONS 

Under  this  heading  of  Recommendations  should  be  listed  the 
solutions  to  the  problems  posed  by  the  Results,  At  the  risk  of 
redundancy  and  speaking  only  for  the  mixtures  in  question  these 
inadequacies  are  repeated  here. 

A  single  equation  of  state  for  all  mixtures  would  be  the 
most  satisfactory  solution  as  it  would  provide  in  one  analytical 
expression  and  its  derivatives  the  measureable  and  derived  prop- 
erties for  the  infinite  number  of  mixtures.   A  solution  of  this 
type  would  require  an  extensive   experimental  program  to  deter- 
mine a  sufficient  number  of  PVT  data  to  permit  the  evaluation 
of  appropriate  coefficients. 

Another  satisfactory  program  would  be  to  make  direct  measure- 
ments of  the  excess  heat  function  using  one  of  the  calorimeters 
mentioned  previously.   This  is  the  only  satisfactory  means  of 
determining  this  property  and  at  high  pressure  the  only  means. 
Then  since  the  vapor  pressure  method  of  determining  G"  is  only 
satisfactory  at  essentially  atmospheric  pressure,  the  excess  en- 
tropy function  could  be  determined  by  the  Integration  procedure 
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discussed  earlier  under  Graphical  Methods,   A  program  of  this  type 
would  not  be  all  inclusive  in  its  results  and  would  necessitate 
taking  measurements  at  sufficiently  close  intervals  of  temperature 
to  permit  reasonable  interpolation  of  tabulated  data  or  accurate 
drafting  of  thermodynamic  diagrams. 
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APPENDIXES 


Appendix  A 

1.  Temperature-composition  diagram,  1  atm. 

2.  Temperature-composition  diagram,  5  atm. 

3.  Temperature-composition  diagram,  10  atm. 

4.  Temperature-composition  diagram,  20  atm. 
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APPENDIX  B 

Barker's  method: 

The  most  widely  noted  method  of  determining  G   from  measured  iso- 
thermal values  of  total  vapor  pressure  and  liquid  composition  is 
that  set  out  by  Barker. 

The  excess  chemical  potential  is  given  by  Eq.  (H-7)    which  can 
be  written 

V-t2   =  RT  ln  N2  "  RT  ln<Pxl£/PN2XN2>  +  (VN2  ~  BN2)(PN2  "  P> 

+  p6n2o2(xo2'2  <66> 

with  a  similar  equation  for  /*-q?» 
It  follows  that 

p  =  7n2pn2  +  To2po2  (6?) 

where  PN  and  Pq  represent  partial  pressures  if  the  liquid  is  an 
•ideal  solution  but  the  vapor  is  a  non  ideal  mixture,   They  are 
given  by 

PN2  =  42PK2  eX?   (VN2  "  BN2"P  "  PN2)/RT  "  P  6n2O2^02'2/RT  <68> 
with  a  similar  expression  for  PQ  and  where 

6n202  =  2BM202  '  BN2  -  B°2  (69) 

If 

°E  =  XN2^N2  +  *02^02  "  42  X02   a  +  ^l2      *02>  +    C^N2-X02)2-'- 

(70) 

then 


8^ 


ln/N2  =  A  XN2  +  B  mN2  +  C  nN2   ..-■'"         .  (71) 

and  similarly  for  In  ffQ    . 

The  values  of  1,  in  and  n  are  expressed  as 

1N2=(X02)  '   102=(x02) 

mN2=-^02)2(1  '  ^42>'     m02  =  (xN2)2d  -  ^o2)  (72) 

nN2=(x02)2  X  -  8xN2  +  12(xN2)2   n02=(xN2)2  1    ~    8X02  +  12(X02)2 

and  A  =  a/RT,  B=b/RT,  C=c/RT. 

As  a  first  approximation  the  solution  is  assumed  to  be  regular, 

such  that  B=C=0.   By  neglecting  corrections  for  the  vapor  phase 

non  ideality,  A  is  set  equal  to  ^ln(2P    /P   +  PN  ), 

P„  and  P   are  then  calculated  for  the  various  experimental 
N2       02 

liquid  compositions  from  Eqs,  (68),   These  are  used  to  calculate 

To  ■  If  11     and  Pcalc  from  Sqs*  (7D  &    (67)  which  in  turn  yield 
the  pressure  residual 

R  =  P     -  P  .  (73) 

expt    calc  y,JI 

and  the  derivatives 

dP/dA  =1   /   p   +  1    T  p"  (7k) 

N2  v   N2  N2     02  U  02  02 

with  similar  expressions  for  dP/dB  and  dP/dC  using  m  and  n  re- 
spectively. 

Changes  in  A,  B  &  C  to  most  nearly  reduce  the  pressure  residuals 
to  zero  are  determined  by  fitting  by  least  squares  to 
(dP/dA)5A  +  (dP/d3)(53  +  (dP/dC)(jC  =  R  (75) 
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or  by  solving  simultaneously  equations  of  the  form 

(5a  2^  (dP/dA)2  +  6b  2_,(dP/dA)(dP/dB)  +  5  C  )_,  (dP/dA)  (dP/dC )  = 

^R(dP/dA)  (76) 

with  similar  expressions  for  )  R(dP/dB)  and  )  R(dP/dC). 


Appendix  C 
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Selected  Data 


1.   Equilibrium  Data 

P      T 

X1 

xv 

H1 

Hv 

S1 

Sv 

Atm    °K 

/£mole  N2 

£mole  N2 

j/rcole 

j/mole 

j/mole  °K 

J /mole  °K 

1   77.36 

100 

100 

823.7 

6793.5 

Q2.84 

165.00 

78 

89.8 

97.3 

722.0 

6790.0 

95.62 

166.23 

80 

64.0 

87.6 

472.0 

6727.0 

98 .  00 

169.18 

82 

44.3 

77.0 

283.0 

6666.0 

98.27 

171.28 

84 

29.6 

63.8 

160.0 

6607.0 

97.88 

173.05 

86 

18.0 

47.6 

70.0 

6547.0 

97.08 

174.00 

88 

8.3 

28.0 

8.0 

6486.0 

95.95 

173.88 

90 

1.0 

3.2 

-29.0 

6426.0 

94.50 

170.73 

90.18 

0 

0 

-32.0 

6406.1 

93.94 

169.63 

5   94.18 

100 

100 

1811.5 

6672.3 

104.21 

155.84 

96 

79.3 

91.0 

I630.O 

6735.0 

108.25 

159.05 

98 

61.3 

80.4 

1472.0 

6800.0 

109.49 

161.33 

100 

45.7 

'69.2 

1337.0 

6865.O 

109.60 

163.09 

102 

33.0 

56.7 

1245.0 

6928.0 

109.28 

164.32 

104 

22.0 

42.4 

1171.0 

6994.0 

108.58 

164.80 

106 

12.3 

26.6 

1108.0 

7053.0 

107.45 

164.58 

108 

3.6 

9.3 

1057.0 

7116.0 

105.78 

162.90 

108.99 

0 

0 

1037.9 

7144.7 

104.60 

160.63 

10  103.95 

100 

100 

2440.2 

6702.3 

110.34 

151.36 

104 

99.3 

99.5 

2433,0 

670^.0 

110.57 

151.53 

106 

79.3 

89.0 

2265.0 

6783.0 

114.27 

155.12 

108 

63.6 

78.7 

2137.0 

6855.0 

115.38 

157.24 

110 

49.6 

68.0 

2028.0 

6926.O 

115.65 

158.88 

112 

37.3 

56.2 

1939.0 

6994.0 

115.49 

159.95 

114 

26c3 

43.4 

1858.0 

7.060.0 

114.89 

160.55 

116 

16.2 

29.1 

1784.0 

7125.0 

113.79 

160.50 

118 

6.9 

13.8 

1721.0 

7186.0 

112.12 

159.30 

119.86 

0 

0 

1681.7 

7245.1 

110.09 

156.50 

20  115.82 

100 

100 

3321.9 

6515.4 

117.96 

145.53 

116 

98o  2 

99.0 

3315.0 

6529.O 

118.48 

146.00 

118 

82.3 

89.1 

3177.0 

6649.0 

121.52 

149.78 

120 

68.8 

79.4 

3072.O 

6752.0 

122.79 

151.98 

122 

56.1 

69.2 

2972.0 

6839.0 

123.09 

153.70 

124 

44.1 

58.0 

2873.0 

6919.0 

122.98 

154.91 

126 

33.0 

46.0 

2785.0 

6993.0 

122.46 

155.62 

128 

22.8 

33.7 

2707.0 

7058.0 

121.53 

155.77 

130 

13.1 

20.7 

2631.0 

7125.0 

120.13 

155.23 

132 

4.1 

7.3 

2543.0 

7190.0 

118.13 

153.77 

133.03 

0 

0 

2542.0 

7208.3 

116.61 

151.69 

2.      Pure    Component   Data 
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p 

T 

H02 

HN2 
j/raole 

S0?  n 

s^2  „ 

Atm 

°K 

j/mole 

j/mole  °k 

j/raole  °K 

1 

65 

-1347.1 

108.4 

76.95 

82.77 

70 

-1115.3 

397.6 

80.38 

87.05 

75 

-861.6 

687.0 

83.88 

91.05 

80 

-594.3 

6486.3 

87.32 

166.02 

82 

-485.3 

6546.8 

88.67 

166.77 

84 

-375.9 

6607.3 

89.99 

167.50 

86 

-265.2 

6667.8 

91.30 

168.21 

88 

-154.0 

6728.O 

92.57 

168.90 

90 

-42.3 

6788.3 

93.83 

169.56 

95 

6940.0 

6938.4 

171.21 

171.20 

100 

7091.1 

7088.1 

172.76 

172.73 

5 

80 

-587.5 

983.1  . 

87.28 

94.68 

• 

85 

-313.5 

1270.8 

90.60 

98.18 

90 

-135.6 

1561.9 

93.78 

101.52 

96 

300.9 

6737.0 

97.40 

156.51 

100 

526.4 

6876.3 

99.70 

157.91 

104 

752.9 

7012.4 

101.92 

159.26 

108 

981.0 

•  7145.8 

104.07 

160.52 

110 

7179.5 

7211.9 

160.94 

161.13 

115 

7349.4 

7374.9 

162.45 

162.56 

120 

7515.3 

7535.7 

163.87 

163.93 

10 

100 

533.9 

2178.0 

99.63 

107.79 

104 

759.8 

6704.5 

101.84 

151.38 

106 

873.2 

6790.2 

102.92 

152.20 

108 

987.1 

6872.9 

103.99 

152.95 

110 

1101.6 

6953.3 

105.04 

153.71 

112 

1216.9 

7031.7 

106.07 

154.41 

114 

1333.1 

7108.5 

107.10 

155.08 

116 

1450.5 

7183.6 

108.12 

155.75 

118 

1569.5 

7257.5 

109.14 

156.37 

120 

7250.7 

7330.7 

156.55 

156.99 

20 

105 

830.3 

2508.0 

102.23 

110.59 

110 

1113.5 

2856.0 

104.86 

113.84 

115 

1400.8 

3249.0 

107.41 

117.32 

118 

1576.5 

6663.6 

108.92 

146.89 

122 

1816.7 

6894.7 

110.92 

148.72 

126 

2066.4 

7096.5 

112.94 

150.35 

130 

2329.7 

7280.5 

115.00 

151.80 

135 

7312.6 

7494.8 

152.46 

153.^0 

140 

7552.3 

7697.1 

154.21 

15^.89 

145 

7770.4 

7890.4 

155.74 

156.23 
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3.   Deviations 


X 

IT 

sM 

SE 

^mole  N? 

j/roole 

j/mole  °K 

J/mole 

0 

0 

0 

0 

10 

36.4 

2.73 

.26 

20 

56.1 

4.14 

.35 

30 

65.7 

5.09 

.38 

ko 

68.6 

5.60 

.37 

50 

66.9 

5.73 

.34 

6o 

62.3 

5.60 

.29 

70 

55.2 

5.09 

.28 

80 

43.9 

4.1*4- 

.24 

90 

26.8 

2.73 

.18 

100 

0 

0 

0 

4-.   Summary  of  equimolar  G 


T 
°K 

GE 
j/mole 

Fig. 33 
notation 

Source  of 
calculation 

Source  of  vapor 
pressure  data 

63.14 

10.88 

A 

7 

7 

65 

10.76 

B 

1 

9 

70 

10.28 

ii 

it 

ii 

77.5 

9.74 

ii 

ti 

ii 

75.16 

9.65 

C 

5 

.8 

77.5 

10.75 

D 

it 

5 

77 

10.0 

E 

2 

2 

83.82 

9.2 

F 

6 

6 

90.0 

9.6 

G 

it 

10 

91.0 

7.3 

it 

ti 

10 

80 

80^5 

H 

This  work 

8 

90 

8.82 

it 

n 

it 

100 

9.75 

n 

it 

it 

110 

9.08 

it 

it 

it 

120 

7.64 

ii 

ii 

ti 

Source  notation  refers  to  "References  for  Excess  Functions" 


